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Abstract. This paper presents the first substantial study of the chemistry of the envelopes around a sample of 18 
low-mass pre- and protostellar objects for which physical properties have previously been derived from radiative 
transfer modeling of their dust continuum emission. Single-dish line observations of 24 transitions of 9 molecular 
species (not counting isotopes) including HCO+, N 2 H+, CS, SO, S0 2 , HCN, HNC, HC 3 N and CN are reported. 
The line intensities are used to constrain the molecular abundances by comparison to Monte Carlo radiative 
transfer modeling of the line strengths. In general the nitrogen-bearing species together with HCO + and CO 
cannot be fitted by a constant fractional abundance when the lowest excitation transitions are included, but 
require radial dependences of their chemistry since the intensity of the lowest excitation lines are systematically 
underestimated in such models. A scenario is suggested in which these species are depleted in a specific region of 
the envelope where the density is high enough that the freeze-out timescale is shorter than the dynamical timescale 
and the temperature low enough that the molecule is not evaporated from the icy grain mantles. This can be 
simulated by a "drop" abundance profile with standard (undepleted) abundances in the inner- and outermost 
regions and a drop in abundance in between where the molecule freezes out. An empirical chemical network is 
constructed on the basis of correlations between the abundances of various species. For example, it is seen that 
the HCO + and CO abundances are linearly correlated, both increasing with decreasing envelope mass. This is 
shown to be the case if the main formation route of HCO + is through reactions between CO and H^", and if the 
CO abundance still is low enough that reactions between Hg~ and N 2 are the main mechanism responsible for the 
removal of H^J". Species such as CS, SO and HCN show no trend with envelope mass. In particular no trend is seen 
between "evolutionary stage" of the objects and the abundances of the main sulfur- or nitrogen-containing species. 
Among the nitrogen-bearing species abundances of CN, HNC and HC3N are found to be closely correlated, which 
can be understood from considerations of the chemical network. The CS/SO abundance ratio is found to correlate 
with the abundances of CN and HC3N, which may reflect a dependence on the atomic carbon abundance. An anti- 
correlation is found between the deuteration of HCO + and HCN, reflecting different temperature dependences 
for gas-phase deuteration mechanisms. The abundances are compared to other protostellar environments. In 
particular it is found that the abundances in the cold outer envelope of the previously studied class protostar 
IRAS 16293-2422 are in good agreement with the average abundances for the presented sample of class objects. 
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1. Introduction 

Understanding the chemistry of protostellar environments 
is important in order to build up a complete and consistent 
picture for the star-formation process. Detailed knowledge 
about the chemistry is required in order to fully under- 
stand the physical processes since it regulates the ion- 
ization balance and the gas temperature through cooling 
of the molecular gas, for example. At the same time the 
chemistry may potentially serve as a valuable tool, both 
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as a time indicator for the protostellar evolution and as 
a diagnostic of the properties of different components of 
young stellar objects. In this paper a molecular inventory 
for the envelopes around a sample of low-mass protostars 
is presented. Compared to previous studies this is the first 
substantial study of the chemical composition of a signif- 
icant sample of low-mass protostars. These objects have 
formed a central protostar, but are still deeply embed- 
ded in their envelope and represent the first stage after 
the collapse of the dark cloud cores. Such objects differ 
from the pre-stellar cores for which chemical stud i es hay e 
previously been performed by, e.g.. iBergin et al ] l|200lh . 
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iTafalla et all l)2002|) and lLee et all l|2003|h in that the cen- 
tral source heats the envelope and dominates the energy 
balance rather than, e.g., the external interstellar radia- 
tion field. 

Other studies of the chemistry of low-mass proto- 
stellar objects include t hose in the Serpens region by 
iHoeerheiide et all l|l999|) and of specific molecules such 
as th e sulfur-bearing speci es and deuterated molecules 
(e.g., iBuckle fc Fuller! Eool iRoherts et ahlEio^. Single 
objects, such as the low-mass protostar IRAS 16293-2422, 



1994; 


van Dishoeck et alJ Il995l 




Ceccarelli et alJ Il998: 


Schoier et alJ l2002t ICazaux et alJ 


120031). This obiect is 



particularly interesting because of its rich spectrum and 
evidence for evaporation of ices in the inner hot regions 
l|Ceccarelli et al.ll2000albT: ISchoier et al.ll2002j) . One of the 
questions that can be addressed with this study is how 
representative the chemistry of IRAS 16293-2422 is com- 
pared to that in other low-mass protostellar objects. 

One of the major steps forward in this line of research 
in recent years has been the observations and analysis 
of the (sub)millimeter continuum emission from the dust 



eras such as SCUBA on the JCMT (e.g 


Chandler & Richer 


200fi iHogerheiide & 


Sandelll 1200(1 


Evans et all l200ll 


Motte AndrJ2()mU 


Jorgensen et all 


2002:IShirlev et alJ 


20021 ISchoier et al.ll2002l) and infrared extinction studies 


Ce.g..lAlves et al.ll200lHHarvev et alJl200l|). Bv fitting the 



radial distributions of the continuum emission and SEDs 
of the objects, the dust component and physical structure 
of the envelopes can be constrained, and, with assump- 
tions about the gas-to-dust ratio and gas-dust temper- 
ature coupling, the physical properties of the gas in the 
envelope can be derived. Such physical models can then be 
used as the basis for determining the molecular excitation 
and for deriving abundances relative to H? by compar- 
ing to molecular line o b servations (e.g..lBergin et all2001t 



Jorge 



nsen et al 



2002; Lee et all 



Line observations (e.g., iSergm et al. zvvi\ 
l2002t ISchoier et all 120021: ITafalla et all 



2003). 



Based on these methods it has become increasingly 
clear that large variations of molecular abundances can 
occur in protostellar environments. Examples are deple- 
tion of molecules at low temperatures due to freeze-out on 
dust grains (e.g..lCaselli et al.ll999tlJorgensen et all2 002: 
Tafall a et al . 2002) and enhancements of molecular species 
in warm regions where ices evaporate ( Ceccarelli et all 
l2000albl : ISchoier et alJ 12002ft or in shoc ked gas associ- 
ated with protostellar outflows or jets jBa^hjUgr^^alJ 
19951 iBachiller fc Perez Gutierrezl fl997t iJorgensen et alJ 



2003 

l.lorgensen et aD l)2002t) (Paper I in the following) es- 
tablished the physical properties of the envelopes around a 
sample of low-mass protostars from ID radiative transfer 
modeling of SCUBA dust continuum maps. The derived 
density and temperature structure and size was used as 
input for modeling CO (sub)millimctcr line emission. It 
was found that the CO lines could be reproduced with 
the physical models assuming constant fractional abun- 



dances with radius. The derived values for the envelopes 
with the most massive envelopes - typically interpreted 
as the "youngest" class protostars - were found to be 
lower than abundances quoted for nearby dark clouds by 
an order of magnitude. In contrast the potentially more 
evolved class I objects were found to have envelopes with 
CO abundances closer to the dark cloud value. It was sug- 
gested that this was related to CO freezing out on dust 
grains at low temperatures and in dense en vironments. 

This paper is a continuation of |Paper"H and the anal- 
ysis of the class YSO, IRAS 16293-2422 presented by 
ISchoier et all l|2002|) . Based mainly on JCMT and Onsala 
20 m observations, abundances for a large number of 
molecules are derived using detailed Monte Carlo radia- 
tive transfer for t he full s et of pre- and protostellar ob- 
jects presented in IPaper J with the envelope parameters 
derived in that paper as input. The combination of low J 
3 mm observations (from the Onsala telescope) and higher 
J lines from the JCMT allows the radial variation of the 
chemistry to be discussed with the low J lines mainly sen- 
sitive to the outer cold part of the envelope and the high 
J lines to the inner dense regions. Similar analyzes for 
H2CO, CH3OH and more complex organic species, which 
are particularly sensitive to the inner most hot core re- 
gion, are presented in separate papers ijMaret et al.l l2003: 
J0rgensen et al. in prep.) 

The paper is laid out as follows: in Sect. [5] the de- 
tails of the observations and reduction are presented. The 
modeling approach is described in Sect. Eland caveats and 
implications for the radial structure described. Relations 
between the different molecular species are discussed in 
Sect.H 

2. Observations 

2.1. Observational details 

The principal data set forming the basis of this work 
was obtained with the James Clerk Maxwell Telescope 
(JCMT) on Mauna Kea, Hawaii 1 where 15 sources were 
observed between February 2001 and February 2003. In 
addition archival data for 3 class I sources - L1551-I5, 
TMC1 and TMC1A - observed previously in a number of 
these settings were used. 

The A3 and B3 receivers at 210-270 GHz and 315- 
370 GHz, respectively, were used with the digital autocor- 
relation spectrometer (DAS) in setups with bandwidths 
ranging from 125 MHz to 500 MHz with resulting reso- 
lutions of 0.1 to 0.6 km s _1 . Each setting was observed 
with on source integration times ranging from 10 to 60 
minutes per mixer reaching a typical RMS (on scale) 
of 0.03 to 0.05 K in 30 minutes. The pointing accuracy for 
the JCMT was found to be a few arcseconds. The calibra- 
tion was checked by comparison to spectral line standards 



The JCMT is operated by the Joint Astronomy Centre in 
Hilo, Hawaii on behalf of the parent organizations: PPARC in 
the United Kingdom, the National Research Council of Canada 
and The Netherlands Organization for Scientific Research 
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and was estimated to be accurate to approximately 20%, 
when comparing data taken in separate runs. For most 
sources beam switching with a chop of 180" was used. 
The only exception was N1333-I4A and -I4B for which 
position switching to an emission-free position at (-120", 
250") was used. 

Further observations at 3 millimeter (85 to 115 GHz) 
were obtained with the Onsala Space Observatory 20 m 
telescope 2 in observing runs in March 2002 and May 2003. 
The entire sample was observed at Onsala in the same 
species, except the two p Ophiuchus sources L1689B and 
VLA1623 which are located too far south. These two 
sources were observed in early April 2003 at 3 mm using 
the Swcdish-ESO Submillimeter Telescope (SEST) 3 at La 
Silla in Chile. Finally CS and C 34 S 2-1 and 3-2 spectra 
were taken for a few sources in November 2001 with the 
IRAM 30 m telescope 4 at Pico Veleta, Spain in the range 
90 to 250 GHz. 

In addition to the observed settings the public JCMT 
archive was searched for useable data and included to con- 
strain the models together with previously published ob- 
servations. All spectra were calibrated at the telescopes 
onto the natural antenn a temperature scale, T * A , using the 
chopper- wheel method l|Kutner fc Ulichlfl98lf) . The spec- 
tra were corrected for the telescope beam and forward 
scattering efficiencies and brought onto the main beam 
brightness scale, T m b, by division with the appropriate 
main beam efficiencies, r\ m \, (or F e e/B c ff in the terminol- 
ogy adopted at the IRAM 30 m telescope). Finally a low 
order polynomial baseline was subtracted for each spec- 
trum. An overview of the observed lines is given in Tabled 

2.2. Resulting spectra 

Spectra of selected molecular transitions are presented in 
Fig. Ill II In order to derive the line intensities, Gaussians 
were fitted to each line. For the few asymmetric lines the 
emission was integrated over ± 2 km s _1 from the sys- 
temic velocity of the given source. The integrated line in- 
tensities are listed in Tables |2"1?)1 In case of non-detection, 
the 2a upper limit is given where u = 1.2%/ Av Sv cr rms 
with Av the expected line width 1 km s _1 for the 
observed sources/molecules), 5v the channel width in the 
given spectral line-setup and <7 rms the rms noise in the 
observed spectra for the specific channel width. The fac- 
tor 1.2 represents the typical 20% calibration uncertainty 



2 The Onsala 20 m telescope is operated by the Swedish 
National Facility for Radio Astronomy, Onsala Space 
Observatory at Chalmers University of Technology. 

3 The SEST is operated by the Swedish National Facility for 
Radio Astronomy on behalf of the Swedish Natural Science 
Research Council and the European Southern Observatory. 

4 The IRAM 30 m telescope is operated by the Institut de 
Radio Astronomie Millimetrique, which is supported by the 
Centre National de Recherche Scientifique (France), the Max 
Planck Gesellschaft (Germany) and the Instituto Geografico 
Nacional (Spain). 



Table 1. Summary of the observed lines. 





Lin© 




_L clcbCUJJ c 




9 1 

z— 1 


07 Q81 n 

y / .yolu 


UoU, OtliO 1 




*3 9 
O Z 


1 A(\ QfiQH 
i4o.yoyu 






^ A 


Z44.yoOO 


TPA/TT TR A A/T 

J U1V1 J_ , 1JTwt_1V1 




7 « 
i — O 


o4Z.ooOU 


J U1V1 1 




9 1 
z — 1 


Q« A 1 9Q 
yu.4izy 


O^O TR A A/T ^P^T 




^ A 


Z41.U101 


J U1V1 1 


1 1 V U 


i— u 


oD. i O4o 


U O U . OtliO 1 




% 9 
o Z 


9fiO 9^^ 
ZDU.ZOOO 


J U1V1 ± 




A ^ 
■TO 


o4o.yyoo 


J V_/1V1 J_ 


I J\ v) 


Q 9 


Z1D. 11ZD 


Tr^A/TT 


IN 2Tl 


1 u 


1 7^7 


UoU, O-CjO X 


n v^in 


A ^ 
4 o 


i^a 

O04.0U00 


TPA/TT 


IT KjXS 


i— u 


O0.04UZ 


UoU, otiiO 1 




*3 9 
o z 


9^Q m 1 s 
zoy.uiio 


TPA/TT 


J_y U1N 


^ 9 


91 7 9*^Sfi 
Zl 1 .ZooD 


TPA/TT 


IT IN U 


i— u 


on AfiQfi 
yu.DDoD 


VJO\ J . OtliO 1 




A ^ 
4 D 


ODZ.D0U0 


7PA/TT 




1 u 


1 1 ^ aqi n 
_i_ 10. 4y iu 


UoU, orjo ± 




3-2 a 


340.2478 


JCMT 


HC 3 N 


10-9 


90.9790 


OSO, SEST 


SO 


2 3 - 1 2 


99.2999 


OSO, SEST 




87 -7 6 


340.7142 


JCMT 


S0 2 


3l,3 — 2o,2 


104.0294 


OSO, SEST 




93,7 — 92,8 


258.9422 


JCMT 



Notes: a Hyperfine splitting observable. 



found by comparing to spectral line standards and obser- 
vations from different nights. 

For most sources the line profiles are quite symmetric 
and can be well-represented by the Gaussians: the main 
exceptions are the strong, optically thick HCN 4-3 and 
CS lines toward especially N1333-I4A and -I4B. The HCN 
and CS lines toward these objects seem to be dominated 
by outflow emission. SO2 is only detected in the low exci- 
tation 3i 3 — 2o,2 line toward N1333-I4A and -I4B, and the 
two objects in p Oph, VLA1623 and L1689B. The higher 
excitation 93,7 — 92,8 line was also observed in a setting 
together with H 13 CN 3-2 but was not detected toward 
any source. Furthermore the high J lines of SO were also 
only detected toward the objects in NGC 1333 and toward 
VLA1623, suggesting a chemical effect. 

Some systematic trends can be seen from the Tables 
and Figures. In general the lines arc signifi cantly weaker 
than those found in IRA S 16293-2422 (|Blake et alJll994l 
Ivan Dishoeck et al1ll995h . Especially for the Class I ob- 
jects in our sample (i.p., L1489 and TMR1) a number of 
usually quite strong lines (e.g., HCN 4-3) were not de- 
tected. The effects of the chemistry are also hinted at by 
comparing, e.g., the source to source variations of the HNC 
4-3 and CN 3-2 spectra. An interesting effect can be seen 
for the deuterium-bearing species: note that the DCO + 
3-2 lines are detected toward the pre-stellar cores but not 
the H 13 CO + 3-2 lines and vice versa for the class I ob- 
jects, clearly indicating a higher degree of deuteration in 
the colder pre-stellar cores. How much of this can be at- 
tributed to excitation and simple mass or distance effects 
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is, however, not clear. IRAS 16293-2422 for example has 
the most massive envelope compared to the other class 
objects and is also located closer than, e.g., the other mas- 
sive sources in the NGC 1333 region. In contrast the class 
I objects by their very definition have the least massive 
envelopes, so the absence of lines toward some of these 
objects may simply reflect lower column densities for the 
observed species for these sources. In order to address this 
in more detail it is necessary to model the ful l line r adia- 
tive transfer as was done in I.Tdrgensen et al and 
ISchdier et all l|2002|) . 



3. Modeling 

3.1. Constant abundances in static models 

To model the ch emical abundances t he same approach 
as in IPaper J and ISchoier et al ] ll2002t was adopted: the 
physical structure for each individual source presented in 
IPaper J was used as input for calculating the molecular 
excitation and the d etailed radiative transfer us i ng th e 
Monte Carlo code of iHoeerheiide fc van der TalJ (|2000|) . 
This co de agrees within the Monte Carlo noise with the 
code of ISchdier et alJ l|2002|) that was used to derive the 
abundances for IRAS 162 93-2422. Both codes were fur- 
thermore benchmarked bv Ivan Zadelhoff et 

ID l|2002l) to- 
gether with other li ne radiative transfer codes. The same 
molecular data as in lSchoier et al.l ( 20j22f) were used. These 
are summarized in the database of Schoier et al. (in prep.). 

A few species, e.g., CN and N2H + , show clear hyper- 
fine splitting of the lines. For the CN 1-0 line the indi- 
vidual hyperfine components can easily be disentangled 
(see Fig. and each of these can be modeled as sepa- 
rate lines with individual excitation rates. In general the 
model fits the individual hyperfine components well, al- 
though in the poorest fits the strongest hyperfine com- 
ponent is overestimated in the modeling. For the CN 3-2 
lines at 340.248 GHz three hyperfine components are over- 
lapping. These transitions are optically thin, however, and 
can therefore be modeled as one line. For N2H+ molecular 
data only exist for the main rotational transitions. Again 
this does not pose a problem if the emission is optically 
thin. For the 1-0 transition this is true for most sources 
as also indicated by the observed ratios of their hyperfine 
components. 

For a given line the resulting sky brightness distribu- 
tion was convolved with the appropriate beam and the 
resulting spectrum compared to the observed one. The 
envelope was assumed to be static and the integrated line 
intensity and line width fitted by varying the abundance 
profile and turbulent line broadening. In the first iteration, 
a constant fractional abundance of each molecule relative 
to H2 was assumed. It is found that most lines are fitted 
well with such a description, except for some of the low 
J 3 mm lines. Abundance jumps, e.g., due to evaporation 
of ice mantles as found for IRAS 16293-2422, are not ex- 
cluded by the present observations. However, the region 
where the ice mantles would evaporate (T > 90 K) is typ- 



Table 9. Inferred abundances for SO2 based on observa- 
tions of the 3i,3 — 2o,2 line from the Onsala 20m and SEST 
telescopes. 



Source Abundance 



L1448-I2 


<6.3x 


10" 


-11 


L1448-C 


<2.2x 


10" 


-10 


N1333-I2 


<2.5x 


10" 


- 10 


N 1333-14 A 


1.8x 


10" 


- 10 


T\T 1 000 TAT) 

IN 1333-14x5 


0.6X 


10 


-10 


L1527 


<2.4x 


10" 


-10 


VLA1623 


1.7x 


10" 


-9 


L483 


<5.3x 


10" 


-11 


L723 


<1.3x 


10" 


-10 


L1157 


<8.4x 


10" 


-11 


CB244 


<1.0x 


10" 


-10 


L1489 


<2.3x 


10" 


-9 


LI 544 


<9.0x 


10" 


-11 


L1689B 


6.3x 


10" 


-10 



Table 13. Inferred abundances for CN using the 1-0 hy- 
perfine transitions and reduced x 2 where applicable. 



Source 


Abundance 


Xrcd 


^lincs 


L1448-I2 


8.3x10" 


10 


0.13 


4 


L1448-C 


2.4x10" 


9 


2.4 


4 


N1333-I2 


2.4x10" 


9 


2.3 


4 


N1333-I4A 


6.1x10" 


10 


1.4 


4 


N1333-I4B 


8.6x10" 


10 


1.5 


4 


L1527 


1.6x10" 


9 


0.93 


4 


VLA1623 


3.5x10" 


9 


2.0 


4 


L483 


3.3x10" 


10 


2.3 


4 


L723 


7.2x10" 


10 


1.8 


4 


L1157 


5.1x10" 


10 


1.1 


4 


CB244 


1.0x10" 


9 


3.7 


4 


L1551 


1.1x10" 


9 


3.6 


4 


L1489 


4.6x10" 


9 




l+(3) 


TMR1 


5.0x10" 


9 




l+(3) 


TMC1 


1.3x10" 


8 


5.1 


4 


TMC1A 


2.1x10" 


9 


3.0 


4 


L1544 


6.8x10" 


10 


1.7 


4 


L1689B 


1.6x10" 


10 


3.9 


4 



ically less than 100 AU (« 0.5-1") for our sources, and 
therefore heavily diluted in the beam. Furthermore, the 
lines presented in this study are predominantly sensitive 
to the material at low to intermediate temperatures in the 
envelope. 

Tables 17^131 list abundances together with the number 
of observed lines and the reduced x 2 for each individual 
source for species for which more than one line was ob- 
served. A summary of the abundances for all molecules 
assuming standard isotope ratios (Table I14|) is given in 
Table El In each of these tables the abundances were 
taken to be constant over the entire extent of the enve- 
lope. 

For a range of the molecules (especially CS and HCO + ) 
the main isotopes are not well suited for determining 
chemical abundances since the lines rapidly become op- 
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Table 2. Line intensities (J Tmb dv) for CS, C 34 S and SO transitions from the JCMT and Onsala 20 m telescope. 



CS CT^S SO 





5-4 


7-6 


2-1 


5-4 


23 - I2 


87 - 7 6 


T -1 A A £71 Tn 

L1448-I2 


0.31 


< 0.18 


0.39 




1.7 


< 0.053 


T -1 A A O /~1 

L1448-C 


2.2 


2.0 


0.47 


0.26 


1.7 


< 0.096 


A T 1 O O O TO 

JN 1333-12 


4.0 


4.9 


IRAM 


0.66 


2.9 


1.7 


AT1 OOO T A A 

N1333-I4A 


7.9 


4.8 


IRAM 


1.25 


7.4 


5.8 


N1333-I4B 


0.3 


4.7 


TT> A TV TO 

IRAM 


0.79 


6.3 


1.8 


L1527 


1.8 


0.45 


TT) A TV JVCl 

IRAM 


< 0.061 


0.42 


< 0.067 


VLA1623 


4.3 


1.6 


SEST 


0.28 


SEST 


0.98 


L483 


4.1 


2.5 


0.38 


0.47 


1.4 


< 0.13 


L723 


4.3 


0.67 


0.18 


0.29 


1.4 


< 0.12 


L1157 


0.93 


< 0.094 


0.20 


< 0.075 


2.1 


< 0.21 


CB244 


1.9 


0.87 


0.16 


< 0.12 


0.87 


< 0.074 


L1489 


0.62 


0.71 


< 0.040 


< 0.098 


0.26 


< 0.083 


TMR1 


0.77 


0.60 




< 0.086 


0.28 


< 0.054 


L1551 


3.2 


3.4 




0.35 


0.59 




TMC1 


0.46 








0.39 




TMC1A 


0.58 








0.12 




L1544 


< 0.14 


< 0.094 


0.22 




0.82 


< 0.13 


L1689B 


< 0.16 


< 0.12 


SEST a 


< 0.028 


SEST a 


< 0.099 



a 3 mm observations from the IRAM 30 m and SEST telscopes given in Table Qj] and |SJ respectively. 

Table 3. Line intensities (/ T MB dv) for the H 13 CO+, N 2 H+, DCO+ and DCN transitions from the JCMT and Onsala 
20 m telescope. 







H la CO + 




N 2 H+ 


DCO+ 


DCN 




1-0 


3-2 


4-3 


1-0 


3-2 


3-2 


L1448-I2 


1.6 


0.89 




10.5 


0.91 


0.13 


L1448-C 


2.0 


1.9 




11.7 


1.9 


0.40 


N1333-I2 


1.8 


2.1 


2.7 


14.2 


0.95 


0.46 


N1333-I4A 


2.3 


1.4 




15.9 


3.2 


0.40 


N1333-I4B 


2.1 


0.57 




13.5 


2.0 


0.31 


L1527 


2.2 


1.1 


0.47 


4.4 


0.70 


0.10 


VLA1623 


SEST a 


4.7 


3.0 


SEST a 


3.6 


0.43 


L483 


1.5 


1.7 


1.2 


13.7 


0.85 


0.25 


L723 


0.61 


0.94 


0.70 


4.8 


0.36 


< 0.064 


L1157 


0.89 


0.61 


0.52 


7.3 


0.64 


< 0.078 


CB244 


0.87 


0.73 


0.43 


6.3 


0.34 


< 0.083 


L1489 


0.96 


0.82 


0.61 6 


0.26 


< 0.092 


< 0.089 


TMR1 


1.1 


0.51 


0.2Q b 


0.37 


< 0.095 


< 0.089 


L1551 




2.4 b 


2.2 b 


15.7 


1.3 




TMC1 




<0.27 6 


<0.12 b 


5.8 


< 0.16 




TMC1A 




<0.27 6 


<0.13 6 


7.1 


< 0.12 




L1544 


0.92 


<0.10 




6.0 


0.56 


0.12 


L1689B 


SEST a 


0.22 




SEST a 


0.55 


< 0.107 



a SEST 3 mm observations of the southern sources given in Table© ''from lHogerheiide et all il997f) . 

tically thick. Moreover the emission from these species in 
the envelope is in some cases hard to disentangle, since 
the line profiles show clear signs of wing emission due to 
outflows and asymmetries attributed to infalling motions 
llGregersen et al.lll997ll200(ilWard-Thompson fc Bucklevl 
200lj). The lines from the weaker isotopes (e.g., C 34 S and 
H 13 CO + ), however, usually do not suffer from these prob- 
lems and were therefore used to constrain the abundances 
where detected. 
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Table 4. Line intensities (/ T MB dv) for the HCN, H 13 CN, HNC and CN transitions from the JCMT and Onsala 20 m 
telescope. 





HCN 




CN 


HNC 






CN 








4-3 


1-0 


3-2 


i— u 


4— o 


1 nn - — 1 )m i 


1 (IT) ( )f 11 O 


1 no i — [ Ini i 


1 noo — ( )m o 


3-2 6 


L1448-I2 


1.4 


0.44 


< 0.13 


o.l 


U.DU 


1.3 


2.7 


0.88 


0.92 


1.2 


T, 1448-1"! 


5.3 


0.82 


0.50 


8.0 


3.4 


1.4 


2.5 


0.87 


1.6 


1.9 


N1 333-T2 


5.3 


0.59 


0.62 


0.9 


O.O 


2.5 


2.8 


1.1 


1.7 


1.5 


N1333-I4A 


6.6 


1.5 


0.65 


11C 
11.0 




1.9 


3.0 


1.2 


1.8 


1.1 


N1 333-T4R 


9.6 


0.96 


0.70 


1 -1 
4.4 


z.U 


1.8 


2.5 


0.88 


1.2 


0.73 


L1527 


1.4 


0.43 


< 0.067 




o.o 


U.4U 


1.4 


2.2 


1.0 


0.81 


1.3 


VLA1623 


2.2 


SEST a 


< 0.153 


D-CjO 1 


Z.O 




SEST" 






1.4 


L483 


3.5 


0.69 


0.52 


O. ^ 


o o 
z.z 


0.63 


1.4 


0.46 


0.99 


2.4 


L723 


1.7 


0.16 


< 0.080 


Z.O 


1 Q 
l.O 


0.28 


0.51 


0.21 


0.31 


1.1 


L1157 


1.1 


0.23 


< 0.087 


3.7 


1.0 


0.57 


0.87 


0.37 


0.35 


0.40 


CB244 


3 9 


0.14 




3.7 


0.74 


n ^4 




n ^1 

U.O-L 


n 4^ 


1.1 


L1489 


0.85 


0.29 


< 0.087 


2.1 


1.0 


< 0.14 


0.63 


< 0.14 


< 0.14 


0.67 


TMR1 


0.63 


0.15 


< 0.15 


2.3 


0.33 


< 0.14 


0.37 


< 0.14 


< 0.14 


0.47 


L1551 




< 0.12 




9.0 


2.97 


1.4 


1.8 


1.2 


1.1 




TMC1 




< 0.094 




3.2 


<0.27 


0.63 


0.72 


0.82 


0.30 




TMC1A 




< 0.092 




3.8 


0.29 


0.55 


0.79 


0.49 


0.49 




LI 544 


< 0.12 


0.70 


< 0.11 


4.0 


0.35 


0.56 


0.90 


0.44 


0.65 


0.21 


L1689B 


< 0.11 


SEST a 


< 0.10 


SEST" 


<0.24 




SEST a 






< 0.092 



Notes: a SEST 3 mm observations of the southern sources given in Table [5] b Hyperfine components blended. 



L1448IRS2 


L144BMM 


N1333I2 0. 


N1333I4A 


N1333I4B 

X 


1.1527 


1483 1 0. 


1723 


LI 157 


CB244 


L1489 1. 


TMRl 1 

fffiiflffl w WI™f 


LI 544 P 


VLA1623 xo.20 

i 


L1689B P. 


C 1B 3-3 



-5 5 10 15 -5 5 10 15 -5 5 10 15 
V [km s- 1 ! 



L1448IRS2 - 


L1448MM - 

^^^^^^^^^^^^^ 


N1333I2 0- 


N1333I4A 










N1333I4B 0- 


14527 0- 


1483 0- 


1723 








L1157 0- 


CB244 1 - 


L1489 L 


TMRl 1 








L1544 P- 


VLA1823 0- 


L1689B P- 


CS 7-6 
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V [km s- 1 ! 



Fig. 1. Spectra of C 18 J = 3-2 (left) and CS J = 7 - 6 (right) from JCMT observations. In this figure, and Fig. 2-gl 
the classes of the individual objects are indicated in the upper right corner of each plot by "0" for the class objects 
(envelope mass > 0.5M©), "I" for the class I objects (envelope mass < 0.5M Q ) and "P" for the pre-stellar cores. 



Table 5. Line intensities ( J Tmb dv) for CS and C 34 S transitions from the IRAM 30 m telescope. 







CS 


C M S 




3-2 


5-4 


2-1 


L1448-C 


2.7 


2.9 


0.45 


N1333-I2 


6.2 


9.4 


1.3 


N1333-I4A 


11.7 


15.1 


1.6 


N1333-I4B 


6.7 


10.3 


0.73 


L1527 


a 


1.8 


0.11 



Notes: a 2a upper limit of C 34 S 3-2 intensity of 0.064 K km s" 1 . 
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Table 6. Line intensities (/ Tmb dv) for the 3 mm observations of the southern sources from the SEST. 









VLA1623 


L1689B 


CN 


— : 

J-022 


— On 1 1 


0.73 


0.13 




lo23 


— Om ■> 


1.2 


0.20 




lo21 


— Om 1 
U()ll 


0.70 


0.16 




lo22 


Ooi2 


0.81 


0.15 


c 34 s 


2-1 




0.41 


0.33 


H 13 CO+ 


1-0 




3.1 


1.2 


H 13 CN 


1-0 




0.98 


0.10 


HC 3 N 


10-9 




0.38 


0.052 


HNC 


1-0 




1.5 


1.1 


N 2 H+ 


1-0 




8.0 


6.0 


SO 


2 3 - 


1 2 


3.1 


2.9 



3 -L1448IRS2 - 
2 

«i(#fL#f^ 


L14481 


M - 


N1333I2 0- 


N1333I4A 


3 N1333I4B - 

2 


L1527 0- 


L4B3 


L723 


3 4,1157 - 

3 

1 


CB244 1 - 


L14S9 1- 


TMR1 1 


3 L1544 P- 

2 

1 

3 fV^ VM A w ^vw»v* 


VLA162 

J 


3 xO.50 - 


L1689B P- 


H 13 C0+ 3-2 



L1448IRS2 0- 


L1448MM - 


N1333I2 


N1333I4A 


N1333I4D 0- 


L1527 0- 


L483 0- 


L723 


11157 0- 


CB244 1 - 


1,1489 1- 


TMR1 1 








1,1544 P- 


VLA1623 xO.50 - 


L1689B P- 

: wv r T^K J a UJ ~jl 1 -Jli-| r rv^ 


DC0+ 3-2 



-5 5 10 15 -5 5 10 15 -5 5 10 15 
V [km s- 1 ! 

1ST) L J 



-5 5 10 15 -5 5 10 15 -5 5 10 15 
V [km s- 1 ! 



Fig. 2. Spectra of H 13 CO+ (left) and DCO+ J = 3-2 (right) from JCMT observations. 
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1,1544 
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L1689B P 


HCN 4-3 
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. 


N1333I4A 


N1333I4B 


L1527 


L483 


. 


L723 
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L1489 1 
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1,1544 P 
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L1689B 


P . 
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Fig. 3. Spectra of HCN (left) and HNC (right) J = 4 - 3 from JCMT observations. 
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Fig. 4. Spectra of CN J = 1 — (left) and J = 3 — 2 (right). The J = 1 — observations are from the Onsala 20 m 
telescope and the SEST (marked with ***), the J = 3 - 2 observations are from the JCMT. 
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Table 7. Inferred abundances for CS and C 34 S and reduced \ 2 where applicable. 



Source 


Abundance 


Xrcd 


^lines 


Lines 6 


CS 


T 1 A A ft TO 






1 
1 


0—4 ^JOlvl 1 ) 




9 c w i n— 9 
Z.O X 1U 


9 R 
Z.O 


9 
Z 


0—4, (— D ^JV^lvll ) 


ATI QQQ TO 


1 Qvl r»— 9 
l.oX 1U 


U.yi 


9 

Z 


0—4 , (— J 


i\ti qoq T/1 A 


4. iXlU 


u.uiy 


z 


K A c 7 fi c nPA/TT^ 
0—4 , 1—0 ^JOMl J 


T\T1 TAR 






O 

z 


5-4 c , 7-6 c (JCMT) 


T 1 ^97 


O.U X 1U 


9 G 

z . y 


yl 
4 


^ a 7 r ^ ipmt r"^n d ~i 

4, i O ^Jvylvll, V . }VJ J 


VT A 1 fi9^ 


o « v i n — 9 


9 f-; 


9 
z 


5-4, 7-6 (JCMT) 




1.5xl0 -9 


i 7 

_L . / 


2 


^-4 7-fi nPMTl 




1.3xlCP 9 


1 'V- 


2 




L1157 


1.9xl0 -10 






5-4 ("JCMT) 7-6 HCMT- nd) 


CB244 


4.0xl0~ 9 


0.22 


2 


5-4, 7-6 (JCMT) 


L1551 


4.1 xlO -10 


0.64 


4 


5-4, 7-6 (JCMT, CSO d ) 


L1489 


2.8xl0~ 9 


2.7 


4 


5-4, 7-6 (JCMT, CSO d ) 
5-4, 7-6 (JCMT, CSO d ) 


TMR1 


l.OxKT 8 


3.2 


4 


TMC1 


6.5xl0~ 9 


7.6 


2 


5-4 (JCMT), 7-6 (CSO d ) 
5-4 (JCMT), 7-6 (CSO d ) 


TMC1A 


4.9X10 -10 


1.9 


2 


L1544 


<2.5xl0~ 10 




1 


5-4 (JCMT) 


L1689B 


<3.0xlO _:LO 




1 


5-4 (JCMT) 


C d4 S 


L1448-I2 


4.1xl0 _ii 




1 


2-1 (OSO) 


L1448-C 


l.lxlO" 10 


1.6 


3 


2-1 (OSO, IRAM), 5-4 (JCMT) 


N1333-I2 


1.4xl0~ 10 


0.50 


2 


2-1 (IRAM), 5-4 (JCMT) 


N1333-I4A 


4.6xl0 -11 


0.10 


2 


2-1 (IRAM), 5-4 (JCMT) 


N1333-I4B 


5.5X10" 11 


4.1 


2 


2-1 (IRAM), 5-4 (JCMT) 


L1527 


1.5xl0 -11 




1+(1) 


2-1 (IRAM), 5-4 (JCMT ; nd) 


VLA1623 


1.8xl0~ 10 


3.0 


2 


2-1 (SEST), 5-4 (JCMT) 


L483 


3.1xl0 -11 


9.2 


2 


2-1 (OSO), 5-4 (JCMT) 


L723 


1.0xl0~ lu 


2.7 


2 


2-1 (OSO), 5-4 (JCMT) 


L1157 


3.7xl0 -11 






2-1 (OSO), 5-4 (JCMT ; nd) 


CB244 


7.2XKT 11 




1+(1) 


2-1 (OSO), 5-4 (JCMT ; nd) 


L1551 


3.7xl0- n 




1 


5-4 (JCMT) 


L1489 


<l.lxl0- 10 




(2) 


2-1 (OSO; nd), 5-4 (JCMT; nd) 


TMR1 


<7.9xl0~ 10 




(1) 


5-4 (JCMT; nd) 


L1544 


3.9xl0 -11 




1 


2-1 (OSO) 


L1689B 


1.2X10 -10 




1 


2-1 (SEST) 



"Number of observed lines; number in parentheses indicate number of lines not detected. ^Observed lines; lines not detected 
indicated by "nd". c Comp lex line profile - intensit y denn ed as line integrated over ±2 km s _1 relative to systemic velocity. 
d CSO line intensities from iMoriartv-Schieven et alJ 1^)95). e The 7-6 line alone corresponds to an abundance of 5xl0~ 10 , the 
5-4 line taken alone to 4xl0 -9 . The latter is in agreement with the results from modeling of the C 34 S lines. 
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Table 8. Inferred abundances for SO and reduced \ 2 where applicable. 



Source 


Abundance 


Xrod 


'"Mines 


Lines 








L1448-I2 


7.0xl0" iu 




l+(l) 


2 3 -l 2 


(OSO), 8 7 -7 6 (JCMT ; nd) 






L1448-C 


1.4xl(T 9 




l+(l) 


2 3 -la 


(OSO), 8 7 -7 6 (JCMT ; nd) 
(OSO), 8 7 -7 6 (JCMT) 6 






N1333-I2 


3.4xl0~ 9 


1.3 


2 


2 3 -l 2 






N1333-I4A 


4.6xl0" 9 


0.66 


2 


2 3 -la 


(OSO), 8 7 -7 6 (JCMT) 6 
(OSO), 8 7 -7 6 (JCMT) 6 






N1333-I4B 


3.0xl0~ 9 


0.82 


2 


2 3 -l 2 






L1527 


1.4xl0~ 10 


1.4 


2+(l) 


2 3 -l 2 


(OSO), 4 3 -3 2 (NRAO a ), 8 7 -7 6 (JCMT; nd) 




VLA1623 


1.2xl0~ 8 


1.6 


3 


23-la, 


6 5 -5 4 (SEST), 8 7 -7 6 (JCMT) 






L483 


2.9xl0~ 10 




1+(1) 


2 3 -la 


(OSO), 8 7 -7 6 (JCMT; nd) 






L723 


2.4xl0~ 9 






23-la 


(OSO), 8 7 -7 6 (JCMT; nd) 






L1157 


1.6xl0~ 9 


0.25 


3+(l) 


2 3 -l 2 


(OSO), 2 2 -li, 43-3 2 (NRAO a ), 8 7 - 


-7 8 (JCMT; 


nd) 


CB244 


9.0xl0~ lu 


3.1 


3+(l) 


2 3 -l 2 


(OSO), 2 2 -li, 43-3 2 (NRAO a ), 8 7 - 


-7e (JCMT; 


nd) 


L1551 


1.9xl0~ 10 


0.95 


2 


2 3 -la 


(OSO), 4 3 -3 2 (NRAO a ) 






L1489 


2.0xl0~ 9 




1+(1) 


2 3 -la 


(OSO), 8 7 -7 6 (JCMT; nd) 






TMR1 


4.1xl0" 9 




1+(1) 


2 3 -la 


(OSO), 8 7 -7 6 (JCMT; nd) 






TMC1A 


2.3 xlO" 10 




1 


2 3 -la 


(OSO) 






TMC1 


4.1xl0~~ 9 


0.70 


2 


2 3 -l 2 


(OSO), 4 3 -3 2 (NRAO a ) 






L1544 


4.8xl0~ 10 




1+(1) 


2 3 -l 2 


(OSO), 8 7 -7 6 (JCMT; nd) 






L1689B 


2.6xl0~ 9 


3.3 


2+(l) 


2 3 -l 2 , 6 5 -5 4 (SEST), 8 7 -7 6 (JCMT; nd) 







a 2 2 -li and 4 3 -3 2 NRAO observations from Buc kle fc Fuller! (|2003l) : see further discussion in text. 6 8 7 -7« line very wide; larger 
uncertainty in derived abundance. 



Table 10. Inferred abundances for H 13 CO + and reduced x 2 where applicable. 



Source 


Abundance 


•2 

Xrcd 


"lines 


Lines 


L1448-I2 


9.9xl0 _i2 ' a 




i 


3-2 (JCMT) 


L1448-C 


1.3xlGT 1Q 


0.21 


4 


1-0 (OSO), 3-2 (JCMT, CSO 6 ), 4-3 (CSO 6 ) 

3-2 (JCMT, CSO 6 ), 4-3 (JCMT, CSO 6 ) 

3-2 (JCMT, CSO 6 ), 4-3 (CSO 6 ) 

3-2 (JCMT, CSO 6 ), 4-3 (CSO 6 ; nd) 

1-0 (OSO), 3-2 (JCMT, CSO 6 ), 4-3 (JCMT, CSO 6 ) 

3-2 (JCMT, CSO 6 ), 4-3 (JCMT) 

1-0 (OSO), 3-2 (JCMT, CSO 6 ), 4-3 (JCMT) 


N1333-I2 


4.7xl0~ 11 ' <1 


0.59 


4 


N1333-I4A 


6.1xl0~ 12 ' a 


1.2 


3 


N1333-I4B 


8.8xl0~ 12,a 


2.7 


2+(l) 


L1527 


8.5xl0~ 12 


0.41 


5 


VLA1623 


2.2xl0" 10 ' a 


2.0 


3 


L483 


2.8xl0 -11 


0.63 


4 


L723 


5.8X10" 11 


0.61 


3 


1-0 (OSO), 3-2, 4-3 (JCMT) 


L1157 


8.4xl0~ 12,a 


0.29 


2 


3-2 (JCMT), 4-3 (JCMT) 


CB244 


7.3X10" 11 


0.53 


3 


1-0 (OSO), 3-2, 4-3 (JCMT) 


L1551 


2.3xl0 -11 


1.3 


2 


3-2, 4-3 (JCMT) 


L1489 


2.5xl0" 10 ' a 




1 


3-2 (JCMT) 


TMR1 


3.9xl0" 10 ' a 




1 


3-2 (JCMT) 


TMC1 


<1.4xHT 10 ' c 




2 


3-2, 4-3 (JCMT C ) 


TMC1A 


<l.lxl0~ u,c 




2 


3-2, 4-3 (JCMT C ) 
3-2 (CSO 6 , JCMT; nd) 
3-2 (CSO 6 , JCMT) 


L1544 


5.6xl0~ 12 ' a 




1+(1) 


L1689B 


1.7xRT 11,a 


0.20 


2 



"Fits i nclud i ng 1-0 line intensities give n in Table fTTI See also discussion in text. 6 CSO measurements from iGreeersen et alJ 
l)l997l EjQOOT) : lGrege£s^rij^_IgYj : r^l)2j | jpjj) - except the 4-3 line towards N1333-I2 (Blake, priv. comm.). c Upper limits on line 
intensities frorn lHogerii^rid'eet^l] 1^)971) . Using their HCO + main isotope 3-2 and 4-3 line intensities, HCO + abundances of 
4.7xl0 -9 and 2.2xl0 _10 are found for TMC1 and TMC1A, respectively. 
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Table 11. Inferred abundances for H 13 CO + including the full set of lines. 



Source 


Abundance 


Xrcd" 


-fobs 6 




L1448-I2 


1.1 X 1U 


ii 


15 


1.58 


0.34 


N1333-I2 


4.0 X 1U 


n 


1.6 


1.77 


0.96 d 


N1333-I4A 


6.5x10" 


12 


6.6 


2.32 


0.38 


N1333-I4B 


9.4x10" 


12 


9.4 


2.10 


0.40 


VLA1623 


2.2x10" 


10 


4.1 


3.09 


1.76 


L1157 


9.2x10" 


12 


8.1 


0.89 


0.17 


L1489 


3.1x10" 


10 


6.0 


0.96 


0.44 


TMR1 


3.9x10" 


10 


9.4 


1.13 


0.40 


L1689B 


1.8x10" 


11 


6.9 


1.24 


0.41 


L1544 


6.7x10" 


12 


16 


0.92 


0.17 



"Reduced x 2 including lines from Table ITuI together with 1-0 lines from the Onsala 20 m telescope or SEST. ^Observed 1— 
line intensity ( f rMB df). c Modeled 1-0 line intensity (J Tmb dv) with abundance from Table tTTJl d See also discussion in 
iJeirgensen et alJ i2003tfl . 



Table 12. Inferred abundances for H 13 CN and reduced \ 2 where applicable. 



Source 


Abundance 


Xrcd 


^lincs 


Lines 


L1448-I2 


2.1xl0 _il 




l 


1-0 


(OSO), 3-2 (JCMT; nd) 


L1448-C 


9.0xl0 -11 


2.5 


2 


1-0 


(OSO), 3-2 (JCMT) 




7.7xl0" n 




1 


3-2 


(JCMT) 


N1333-I2 


2.5xl0" n 


7.1 


2 


1-0 


(OSO), 3-2 (JCMT) 




2.9xl0" n 




1 


3-2 


(JCMT) 


N1333-I4A 


5.6xl0 -12 


18 


2 


1-0 


(OSO), 3-2 (JCMT) 




5.1xl0" 12 




1 


3-2 


(JCMT) 


N1333-I4B 


3.2xl0 -11 


1.9 


2 


1-0 


(OSO), 3-2 (JCMT) 




2.8xl0" n 




1 


3-2 


(JCMT) 


L1527 


3.2xl0" n 




1 


1-0 


(OSO), 3-2 (JCMT; nd) 


VLA1623 


2.5xl0 -10 




1+(1) 


1-0 


(SEST), 3-2 JCMT; nd) 


L483 


2.4xlO -11 


1.0 


2 


1-0 


(OSO), 3-2 (JCMT) 




2.9xlO -11 




1 


3-2 


(JCMT) 


L723 


3.6X10" 11 




1 


1-0 


(OSO), 3-2 (JCMT; nd) 


L1157 


2.0x10"" 




1 


1-0 


(OSO), 3-2 (JCMT; nd) 


CB244 


3.1xl0 -11 




1 


1-0 


(OSO), 3-2 (JCMT; nd) 


L1489 


<1.7xl0- 10 




(2) 


1-0 


(OSO; nd), 3-2 (JCMT; nd) 


TMR1 


3.7xl0" 10 




1+(1) 


1-0 


(OSO), 3-2 (JCMT; nd) 


L1544 


7.1xl0" n 




1+(1) 


1-0 


(OSO), 3-2 (JCMT; nd) 


L1689B 


1.3xl0 -11 




1+(1) 


1-0 


(SEST), 3-2 (JCMT ; nd) 



See description in Table |7| 
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3.2. Shortcomings of the models; drop abundance 
profiles 



For a number of species the constant fractional abundance 
model gives poor results (x 2 > 5) when fitting both the 
lowest rotational lines from the Onsala 20 m and higher 
excitation lines from the JCMT. A similar trend was seen 
in modeling of the CO isotopic species in Paper I, where 
the 1-0 lines were typically underestimated in models fit- 
ted to the 2-1 and 3-2 lines. This trend is particularly 
pronounced for H 13 CO + and the nitrogen-bearing species 



(HCN, H 13 CN, CN and HNC), whereas the low J lines of 
CS, for example, can be fitted well by a constant abun- 
dance. This may be due to the critical density of the ob- 
served transitions which should be compared to the typical 
freeze-out and desorption timescales for the given densi- 
ties and temperatures. Fig. [S] shows the density for two 
objects, N1333-I2 and TMR1, as function of temperature 
(i.e., depth) compared to the critical densities of vari- 
ous t ransitions of CS, CO, HCO + and HCN fe.g.. l.Iansenl 
1995) for the same temperatures. Since the critical densi- 
ties of the CS/C 34 S 2-1 lines are higher than those of the 
HCO + and CO 1-0 lines, CS is less sensitive to the outer 
region of the envelope where depletion and contribution 
from the surrounding cloud may be important. This may 
explain why these transitions can be modeled in the con- 
stant abundance framework. The observed 4-3 transitions 
of, i.p., HCN and HNC have the highest critical densities 
and these lines therefore probe the innermost part of the 
envelope. 

In the outer regions of the envelope the depletion 
timescale for CO is comparable to the lifetime of the pro- 
tostars (~ 10 4 — 10 5 years) at the temperatures where the 
molecule can freeze-out. This could explain the failure of 
the constant abundance models in describing the lowest 
J lines for CO (and thereby also H CO + ; see discussion 
in Sect. 14.31) : in prestellar cores (e.g ICaselli et all 1 19991 
iTafalla e t al. 2002) a trend is seen of decreasing CO abun- 
dances with increasing density toward the center. Since 
the temperature in the bulk of the material in these ob- 
jects is low enough for CO to be frozen out, the expla- 
nation for the radial dependence is a difference in density 
and thus the freeze-out timescale. Therefore the time for 
CO to freeze-out in the outermost regions may simply be 
too long to result in appreciable amounts of depletion. For 
the protostellar cores the difference is the heating by the 
central source, which induces a temperature gradient to- 
ward the center. CO is therefore expected to be frozen out 
in a small region, where the density is high enough that 
the freeze-out timescale is short, yet the temperature low 
enough that CO is not returned to the gas-phase. 

A simple way of testing this can be performed by intro- 
ducing a "drop" chemical structure as illustrated in Fig. El 
with a constant undepleted CO abundance Xq in the parts 
of the envelope with densities lower than 3 x 10 4 cm~ 3 or 
temperatures higher than 30 K. A lower CO evaporation 
tem perature of ~ 20 K is ruled out by the 3-2 line intensi- 
ties l|Paner tf l. The undepleted abundance, Xq, is taken to 
be the same in the inner and outer regions of the envelope 
to avoid adding another free parameter. The abundance 
in the region with temperatures lower than 30 K and den- 
sities higher than 3 x 10 4 cm -3 , X&, can then be adjusted 
to fit the observations. 

Fig.0shows a comparison for L723 between the models 
with a constant fractional abundance from Paper I and a 
model with two abundance jumps described above. The 
latter model has two free parameters (besides the Doppler 
broadening, which does not alter the results), Xq and Xd. 
In the constant abundance model, the C 18 abundance 
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Fig. 5. Density as function of temperature for the en- 
velopes around TMR1 and N1333-I2 (solid line) compared 
to the critical densities of the observed transitions of CS, 
CO, HCO + and HCN. The critical densities are indicated 
in order of increasing excitation by the dashed-dotted, dot- 
ted and dashed lines, respectively, i.e., showing the 2-1, 
5-4 and 7-6 transitions for CS, the 1-0, 2-1 and 3-2 tran- 
sitions for CO, and the 1-0, 3-2 and 4-3 transitions for 
HCO+ and HCN. 
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Fig. 6. Simulated abundance profile in "drop" models. 



Fig. 7. Fitted C 18 line-profiles for L723. Upper pan- 
els: constant fractional abundance of 3.9xl0 -9 from 
Paper I. Lower panels: drop model with abundances of 
Xo=2.0xl0 -7 for the undepleted material (with either 
temperatures higher than 30 K or densities lower than 
3xl0 -4 ) and an abundance of X^—l.^y. 10~ 8 for material 
with CO frozen out. 



is 3.9xl0 -8 , while in the "drop" model, the undepleted 
abundance Xq is 2xl0~ 7 and the depleted abundance Xn 
is 2xl0 -8 . Similar fits to the C 18 abundances of one of 
the class I objects, L1489, provide equally good results - 
again allowing the 1-0 lines to be fitted together with the 
2-1 and 3-2 lines. The fitted abundances in the case of 
L1489 are X of 5xl0~ 7 and X u of 5xl0~ 8 . 

The fact that the 1-0, 2-1 and 3-2 lines can all be 
fitted in the drop models is not unexpected since an ex- 
tra free para meter is introduced compared to the results 
presented in IPaper 4 , which is used to fit only one ex- 
tra line. Still, it should be emphasized that the chemical 
structure in the drop models has its foundation in results 
from the pre-stellar cores and is thus not completely ar- 
bitrary. As expected, the constan t fraction al abundances 
found for both L723 and L1489 in IPaper"! are a weighted 
average of Xu and X from the drop models. While the 
constant abundances were significantly different for L723 
and L1489 (1.9xl0~ 5 and l.OxlO" 4 , respectively), those 
in the drop models are more similar: the factor 2.5 dif- 
ference in derived abundances can be explained through 
the uncertainties and approximations in the physical and 
chemical description. 

Besides these limitations, the constant fractional abun- 
dances characterize the overall envelope chemistry as a 
good first-order approximation. The more general trends 
in abundance variation from source-to-source can thereby 
be used for a statistical comparison with the caveat that 
the selected transitions may be probing different tempera- 
ture and density regimes in which the chemistry may vary. 
It is important to note that none of the abundances are 
correlated with the distances to the sources or the slopes 
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of their density profiles, indicating that the uncertainties 
in these parameters do not introduce signficant systematic 
errors. 

VLA1623 shows high abundances of most molecular 
species compared to the average class objects. As men- 
tioned in Paper I the envelope model of this particular 
object is highly uncertain since it is located in a dense 
ridge of material and molecular tracers with low critical 
densities, in particular the CO lines and the low J 3 mm 
transitions of the other species, may be sensitive to this 
component rather than the envelope itself. 

TMC1A stands out among the remainder of the 
class I obj ects with significantly low er abundances in all 
molecules. iHogerheiide et aP |l998) likewise found that 
the envelope mass estimated through 1.1 mm continuum 
observations was a factor 5 higher than the mass estimated 
on the basis of 13 CO, C 18 and HCO + measurements. 
One possibility is that TMC1A does have a more massive 
envelope and thus lower ab undances due to depletion such 
as seen for CO in Paper J. Alternatives could be that the 
density in the envelope of this object has been overesti- 
mated from the models of the dust continuum emission or 
that the molecular line emission is tracing material not di- 
rectly associated with the bulk material in the protostellar 
envelope. 

This could be a general problem for more sources: are 
there systematic errors of the envelope dust mass lead- 
ing to false trends in abundances? A systematic overes- 
timate of the mass (i.e., density) for the class objects 
would lead to systematically lower abunda nces, sim ilar 
to the depletion effects observed for CO in IPaner 1 On 
the other hand a change in abundance as seen, e.g., for 
CO, would require that the density scale for the class 
objects is off by approximately an order of magnitude, 
and the submillimeter dust emission and molecular lines 
would have to trace quite unrelated components. This 
is contradicted by the relative success of the models in 
simultaneously explaining observations of both line and 
continuum emissio n from single-dish telescopes ijPaper J . 
ISchoier et aj] 2002. this paper) and higher r esolution inter- 
ferometers ij.lorgensen et al.ll2003bt ISchoier et al.ll2flo3) . 

3.3. Effect of velocity field 

Most observed lines are simple Gaussians with typical 
widths of 1 km s _1 (FWHM). Still, for some molecules sig- 
nificant variations are found between the widths for differ- 
ent rotational transitions and thereby the broadening due 
to systematic and/or turbulent motions required to model 
the exact line profiles. This indicates either systematic in- 
fall in the envelope as expected from the line profiles of 
some of the optically thick species or a variation of the 
turbulent velocity field with radius. 

The problem with the current models is that the 
power-law density profile adopted in lPaperH does not give 
direct information about the velocity field, as would be ob- 
tained by fitting a specific collapse model like the inside- 



out collapse model bv IShul l)l977j) . A velocity field can, 
however, still be associated with the derived density dis- 
tribution, using the mass continuity equation. This equa- 
tion: 



becomes for a spherical symmetric envelope: 
d P , 1 9 ( 2 \ n 



(1) 



(2) 



where v r is the radial velocity. Assuming p cx r -p , a power- 
law velocity distribution v cx r~ q and a static envelope 
density distribution (|£ = 0) results in: 



r A at 



r 2 p q = const. ^ q = 2 - p 



(3) 



(4) 



So for a given power-law density distribution, n(H2) = 
"na{ r / r o)~ v , it is possible to introduce a correspond- 
ing power-law distribution for the velocity field, v — 
vo(r/ro)~ q . Here a characteristic infall velocity vq is in- 
troduced as an additional free parameter, which can be 
fitted by comparison of the line profiles to the turbulent 
linewidth. 

In Fig. [S] such a comparison is shown for the C 34 S ob- 
serv ations for the "typical" class object, N1333-I2 (see 
also lJorgensen et~ . 2003b). The observed line widths are 
seen to constrain the velocity field in terms of the com- 
bination of turbulent broadening and magnitude of the 
systematic velocity field. For a parameterization of the 
velocity field, an estimate of the mass accretion rate M 
can be derived from: 



M = An r [i mn no f o 



(5) 



where [i is the mean molecular weight, 2.33. For N1333-I2 
the upper limit on vq (at the inner radius, ro = 23.4 AU) 
of 2.5 km s , i.e., assuming no turbulent broadening, 
translates to a mass accretion rate of 3xl0 -5 M© yr . 
This agrees with typical mass accretion rates inferred for 
the youngest protostars (e.g.. IShulll977t iBontemos et alJ 
1992 iDi Francesco et all 12001^ . The advantage of using 
the optically thin species to constrain the velocity field is 
that they do not suffer from confusion with, e.g., outflows, 
but only pick-up the bulk envelope material as illustrated 
by high angular resolu tion interferometer studies (e.g., 
I.Torgensen et al.ll2003bl) . On the other hand, complemen- 
tary information about the velocity field is obtained from 
the detailed line-profiles of the optically thick, strongly 
self-absorbed lines, e.g., the relative strength of red- and 
blue pea ks and the de pth of the self-absorpti on feature 
(see, e.g.. lEvansl l|l999h and lMvers et alJ l|2000l) for recent 
reviews of this topic). 
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Fig. 8. Modeling of the velocity field in the envelope 
around N1333-I2: in the upper panel C 34 S model lines are 
compared with observati ons for a constant broadening of 
0.8 km s _1 , as in |Paner"3 . In the lower panel, a model with 
no turbulent broadening, but a power-law velocity field 
with v = 2.5 km s _1 at the inner radius, r = 23.4 AU, is 
adopted. In both plots a constant abundance of 1.4x 10 -1 ° 
was assumed. 



Fig. 9. Dependence of the C 34 S abundance on velocity 
field for an infalling envelope around N1333-I2. The grey 
region indicates la confidence and the almost vertical lines 
indicate the 2a, 3a and 4cr confidence levels. It is seen that 
the derived abundances do not depend on the assumed 
velocity field for this optically thin species. 

Table 14. Adopted isotope ratios. 



Fig. illustrates the important point that the derived 
abundances do not depend critically on the adopted ve- 
locity field for optically thin species like C 34 S, illustrat- 
ing that the static envelope structure provides an ade- 
quate description to derive their overall chemical proper- 
ties. Th is is in agreement with the conclusion reached in 
IPaper J . Note that the confidence levels on the derived 
abundance in Fig.[!|]only correspond to the calibration er- 
ror. Systematic errors due to uncertainties in the adopted 
model, collisional data etc. are not taken into account, so 
the abundances derived may still be subject to uncertain- 
ties not apparent from this figure. 

In summary, although systematic velocities probably 
exist in all envelopes besides the turbulent motions, the 
derived abundances do not depend critically on the de- 
tailed treatment of the velocity field as long as predomi- 
nantly optically thin lines are considered, as in this paper. 
We will therefore for the remainder of this paper stick with 
the assumption of a non-infalling envelope with a constant 
turbulent broadening reproducing the approximate width 
of the lines. 

4. Discussion 

4.1. General trends and empirical correlations 

A direct comparison of the derived abundances for the 
main isotopes can be found in Table ITSl Where possible 
the abundances are calculated using the optically thin iso- 
topic species with the standard isotopic ratios listed in 
Table El 

In general the derived abundances vary by one to two 
orders of magnitude over the entire sample. Following the 



Isotope ratio Value Reference 
rl G/ m C 70 Wil son fc Rood (19941 

16 0/ 18 540 Wilson fc Rood (19941 

18 0/ 17 3.6 Penzias (1981).Paper I" 

32 S/ 34 S 22 Chin et al. (1996) 

a The 18 0/ 17 ratio is not used for the abundances d erived in 
this paper, but was used for the CO abundances in IPaper J 
and is therefore included here for completeness. 

trend seen in Paper I of increasing abundances with de- 
creasing envelope masses, the objects are accordingly sep- 
arated into groups with envelope masses (M>iok) higher 
or lower than 0.5 M©, roughly corresponding to class and 
class I objects, respectively. This definition only moves the 
two borderline class 0/1 objects L1551 and CB244 from 
class I to class and vice versa compared to the source 
list given in Table 1 of Paper I. 

On average the class objects have lower abundances 
than the class I objects for most species (see Fig, imp . The 
most pronounced effect is seen for CO, HCO + and CN 
where the average abundances differ by up to an order of 
magnitude, whereas especially SO and HCN have close to 
constant abundances with envelope mass, albeit with large 
scatter around the mean. As discussed in the following 
sections the variations of abundances with mass are not 
identical, however, which indicates chemical effects regu- 
lating the relative abundances for the different molecular 
species. In order to quantify this more rigorously and in 
an unbiased way, the Pearson correlation coefficients were 
calculated for each set of abundances and are listed in 
Table [TBI The Pearson correlation coefficient is a measure 
of how well a (x, y) data set is fitted by a linear correlation 
compared to the spread of (x, y) points. Values of ±1 in- 
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Table 15. Overview of derived abundances for main isotopic and deuterated species. 



Source 


CO 


CS 


SO 


HCO+ 


DCO+ 


N 2 H 4 


HCN 


DCN 


HNC 


CN 


HC 3 N 




xl(T 


' 5 xl0~ 9 


xlO" 9 


xlO" 9 


xlO" 11 


xl0~ 


9 xlO" 9 


xlO -11 


xlO- 10 


xlO" 10 


xlO" 10 










Class 


(envelope 


mass > 


O.5M ) 










L1448-I2 


0.61 


0.90 


0.70 


0.69 


0.79 


>1.0 


8.0 


0.33 


0.35 


1.8 


1.9 


L1448-C 


3.7 


2.4 


1.4 


9.1 


9.8 


3.9 


5.4 


4.9 


13 


20 


12 


N1333-I2 


2.4 


3.1 


3.4 


3.3 


1.6 


5.0 


2.0 


2.1 


1.8 


3.0 


4.3 


N1333-I4A 


0.79 


1.0 


4.6 


0.43 


1.2 


>1.0 


0.36 


0.31 


0.28 


0.37 


0.72 


N1333-I4B 


1.3 


1.2 


3.0 


0.62 


2.5 


3.2 


2.0 


1.0 


1.4 


1.4 


1.1 


L1527 


3.9 


0.33 


0.14 


0.60 


2.9 


0.25 


1.2 


1.5 


3.2 


24 


8.9 


VLA1623 


16 


4.0 


12 


15 


17 


>3.0 


6.6 


4.7 


10 


8.9 


3.8 


L483 


1.4 


0.68 


0.29 


2.0 


1.1 


0.75 


2.0 


0.94 


3.9 


8.3 


1.8 


Li 16 


1 . 9 


1.1 


O A 

1A 


4. 1 


1.0 


1 9 

1.3 


1.0 


<0.91 


r 1 


o c 

O.U 


O 7 
1. i 


Til CC7 


n £9 
U.Oz 


U.ol 


1.6 


u.oy 


one; 

u.yo 


v i n 
>1.U 


U.UOO 


^-n 9Q 


U.ol 


U.OO 


l.O 


T 1 rri Tr 

-Llool-lo 


q n 
O.U 


A O 1 

U.ol 


0.19 


1.0 


1 r 


3.1 






U. 11 




O 1 

z. 1 


116293-2422 


3.3 


3.0 


4.4 


1.4 


1.3 


0.14 d 


1.1 


1.3 


0.69 


0.80 


1.5 










Class I 


(envelope 


mass < 


O.5M ) 










L1489 


10 


2.8 


2.0 


18 


< 2.3 


0.15 


0.65 


<6.4 


13 


22 


8.7 


in i 


20 


10 


4.1 


27 


< 5.0 


0.35 


1.6 


<15 


8.1 


47 


35 


r T\\C*'\ A 
1 lvlL>li\ 


Z . o 


u.4y 


U.Zo 


n 99 
u.zz 


< 0.65 


3.9 






U.oo 




1UU 




zu 


. 


A 1 


A 7 


< 8.5 


>1.0 










1 




9 7 
O. i 


1.0 


n on 


0.1 


2.1 


2.0 


4.y 


i n 

<i.y 


7 O 


on 
zU 


u.y 












Pre-stellar 












LI 544 


0.49 


0.86 


0.48 


0.39 


2.1 


5.0 


<0.35 


0.77 


12 


4.8 


16 


L1689B 


2.4 


2.6 


2.5 


1.2 


2.4 


0.43 


<0.38 


<2.1 


<4.6 


<2.3 


0.36 


Averages: 
























'class 0' a 


2.1 


1.5 


2.0 


2.2 


2.3 


2.5 


1.3 


1.4 


2.8 


6.9 


3.5 


'class I' 


11 


4.3 


2.3 


11 


2.1 6 


1.6 


2.1 




7.4 


30 


31 


Pre-stellar 


1.4 


1.8 


1.5 


0.80 


2.3 


2.7 


< 0.36 


0.77 


12 


4.8 


8.2 



"excluding VLA1623. Wv detected for CB244. fechoier et al.l (120021) . d This paper. 



Table 16. Pearson correlation coefficients for the abundances for all objects. 





Mass 


CO 


HCO+ 


CS 


SO 


HCN 


HNC 


CN 


HC 3 N 


Mass 




-0.74 


-0.51 


-0.46 


-0.18 


-0.11 


-0.40 


-0.71 


-0.71 


CO 


-0.74 




0.79 


0.69 


0.35 


0.46 


0.52 


0.69 


0.59 


HCO+ 


-0.51 


0.79 




0.80 


0.48 


0.44 


0.70 


0.69 


0.48 


CS 


-0.46 


0.69 


0.80 




0.79 


0.29 


0.48 


0.31 


0.39 


SO 


-0.18 


0.35 


0.48 


0.79 




-0.05 


0.14 


-0.27 


-0.03 


HCN 


-0.11 


0.46 


0.44 


0.29 


-0.05 




0.63 


0.55 


0.45 


HNC 


-0.40 


0.52 


0.70 


0.48 


0.14 


0.63 




0.86 


0.72 


CN 


-0.71 


0.69 


0.69 


0.31 


-0.27 


0.55 


0.86 




0.83 


HC 3 N 


-0.71 


0.59 


0.48 


0.39 


-0.03 


0.45 


0.72 


0.83 





dicate good correlations (with positive or negative slopes) 
whereas a value of indicates no correlation. 

As can be seen from Table El significant differences 
exist between the various sets of abundances. Setting an 
(arbitrary) cut of > |0.7| to indicate good correlation, 
the results suggest that the molecular species are related 
as indicated in Fig. ^2 Individual results are shown in 
Fig. I12I21A The abundances of groups of species, e.g., the 
nitrogen- or sulfur-bearing species, are closely related as 
expected from naive chemical considerations. HCN is the 
only molecule whose abundance does not directly corre- 
late with that of any other molecule at this level. The 
closest correlation is found with its isomer HNC (corre- 



lation coefficient of 0.63). The best correlation between 
abundance and mass is found for CO followed by CN and 
HC3N. Naturally the correlations seen in this comparison 
may indicate the structure of the general chemical net- 
work rather than direct relations between the individual 
molecules: for example the ranking of correlations for SO 
is as follows: CS (0.79), HCO+ (0.48) and CO (0.35). As 
indicated in Fig. ^] this is exactly the decline in correla- 
tion coefficients one would expect with the relations be- 
tween these species on the pair-by-pair comparison basis 
adopted when constructing Fig. Such "connectivity" 
could also be the cause for the relation between HNC, CN 
and HC3N - the correlation between HNC and HC3N may 
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Fig. 10. Comparison between a verage abundances for class and I objects and pre-stellar cores (this paper), 
IRAS 16293-2422 out er envelope llSchoier et al J 120(121) . average abundances for W3(IRS4), W3Q RS5) and W3(H20 ) 
(all high-mass YSOs: iHelmich fe van Dishoecklll997|) and abundances in the dark cloud L13 4N foickens et al-lfeoOCj) . 
Note that the L1 34N abundances have been resca led assuming a CO abundance of 2.7xl0 -4 l|Lacv et al.ll994[) . as was 
also assumed bv IHelmich fe van Dishoeckl l|l997|) for the high-mass YSOs. The L134N abundances thereby become: 
[CS] = 2.7xl0" 9 , [SO] = 1.5X10" 8 , [HCO+] = 2.1xl0" 8 , [HCN] = 2.0xlCT 8 , [HNC] = 7.0xlCT 8 , [CN] = 1.3xl0- 9 , 
and [HC 3 N] = 1.2xlQ- 9 . 



in fact just reflect that both these molecules are related 
to CN. The rather low number statistics imply that care 
should be taken not to overinterpret the absolute values 
of the correlation coefficients, but as a first step they give 
valuable hints. To fully understand the underlying chem- 
istry, a more in-depth consideration on a species by species 
basis is required, as discussed in the following sections. 

4.2. CS and SO 

As can be seen from Fig. El abundances of the sulfur- 
bearing species, CS and SO, are close to constant with 



envelope mass, contrasting the picture for CO l|Paper J) . 
CS has often been used to c onstrain the density scal es in 
protostellar envelopes (e.g., Ivan der Tak et alT EoOO^l as- 
suming the chemistry to be homogeneous throughout the 
envelope. SO also does not show any significant trends 
with envelope mass, but has a larger scatter. 

CS and SO are indeed found to show relatively 
small abundance variat ions in pure gas-pha se mod- 
els of pre-stellar cores l|Bergin fc Langerl Il997|) and in 
cold exterior regi ons of high-ma ss protostellar envelopes 
llDotv et alJl2002|) . As a rgued bv lBergin fc Langerl l(l997l) 
and iBergin et alJ (|200ll) . however, sulfur-bearing species 
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Fig. 11. Relations between different molecules as judged 
from the Pearson correlation coefficients. The dashed line 
between HCN and HNC indicates the strongest correlation 
for HCN with any of the other molecules studied. The 
correlation coefficient for this relation is, however, lower 
than the cut of 0.7 adopted for good correlations. 

such as CS and SO should suffer from depletion at 
densities and temperatures characteristic for these re- 
gions. Maps of CS toward pre -stellar cores l|Tafalla et alJ 
20021 iDi Francesco et al.| 2002l) and molecular clouds (e.g., 
IC 5146; Be rgin et al.ll200lj) combined with models of the 
abundances suggest that this molecule does indeed freeze 
out toward the inner colder and denser parts. Typical 
abundances in such environments range from «lx 1CP 10 
to a few x 10~ 9 between the inner (low abundance) and 
outer (high abundance) regions. This agrees well with the 
average abundances found for the protostellar envelopes 
analyzed in this paper, which have a central source of heat- 
ing. Our CS abundances are also simil ar to those inferred 
for a s ample of high-mass protostars bv lvan der Tak et alJ 
( 2000) using a similar analysis. 

An important conclusion regarding the derived CS 
abundances concerns the impact of outflow processing 
of the gas in the envelopes: CS and SO are seen to 
be gr eatly enhanced in shocked gas in protostellar out- 
flows llBachiller fe Perez Gutierrezlll997t1jOrgensen et alJ 
2003ah . The small source-to-source variation in the de- 
rived CS abundances, however, illustrates that although 
increased CS abundances may be present in small parts 
of the envelopes, the bulk of the emission originates in 
parts of the envelope unaffecte d by such processes. The 
same conclusion was reached by lj0rgensen et alJ l|2003b|) 
from millimeter interferometer observations of the C 34 S 
2-1 line emission toward NGC 1333-IRAS2. 

It has been suggested that comparison between sulfur- 
bearing species like SO and CS can be used as chem- 
ical probes o f the evolutionary stages in star-forming 
regions (e .g.. iRuffle et all Il999l) - both when consider- 
ing high- llCharnlevlll997t lHatchell et alJll99^ and low- 
mass stars l|Buckle fc Fullerll2003l) . The time-dependence 
of the sulfur-chemistry network is initiated when signif- 
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Fig. 12. Abundances of CS from optically thin C 34 S iso- 
topic lines (where detected) and CS lines (upper panel) 
and of SO (lower panel) vs. mass. In this figure and in 
following figures in this paper, the class objects are in- 
dicated by the class I objects by "0" and the pre- 
stellar cores by The class objects VLA1623 and 
IRAS 16293-2422 have been singled out by "□" and "*" , 
respectively. 

icant amounts of H2S are released in the gas-phase by 
evaporation of grain-mantles. This is followed by forma- 
tion of SO and SO2 (through reactions with H and H30 + 
forming S and H3S + and subsequently through reactions 
with OH and O2). At later times most of the sulfur is 
incorporated int o CS, H?CS and PCS. 

In particular iBuckle fc Fulle"rl l)2003h estimated abun- 
dances of sulfur-bearing species from SO, SO2 and H2S 
line observations toward a sample of class and I objects 
assuming LTE and a constant CO/H2 abundance ratio. 
They found that their class I low-mass YSOs had lower 
abundances of SO and H2S than class objects, suggesting 
that this was a result of their later chemical evolutionary 
stage. For the sources in this paper it is seen, however, 
that there is no significant difference in SO abun dances 
between class and I objects. Van der Tak et al. (2003) 
surveyed a range of different sulfur-bearing species toward 
a sample of high-mass YSOs and likewise found no system- 
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atic trends between known indicators of the evolutionary 
stage and the abundances of the sulfur molecules. 

Th ere is a slight overlap between the objects stud- 
ied by iBuckle fc Fuller! (2003) and thos e treated in this 
paper. For these overlapping objects the IBuckle fc Fuller! 
data have been included in our analysis as indicated in 
Table |SJ and it is found that our models can explain their 
observations. A likely explanation for the different findings 
is therefore the CO depletion found for sources with the 
more massive envelopes (Paper I and Sect. l3~2l in this pa- 
per). In fact abundances calculated assuming a [CO /Ha] 
abundance of lxl0~ 4 leads to overestimated abundances 
for objects in which C O is depl eted, i.e., those with the 
most massive env elopes llPaner Jl . The ab undances for the 
class objects in IBuckle fc Fuller! l)2003h could therefore 
be overestimated and their evolutionary trend an artifact 
of this assumption. Fig. ED compares the relation between 
CS and SO abundances relative to the density scale set in 
Paper I and to a CO abundance of 10~ 4 . Fixing the CO 
abundance increases the average SO and CS abundances 
for the class objects - to almost an order of magnitude 
higher than t hose for the clas s I obj ects. This in fact re- 
sembles what lBuckle fc Fuller! l|2003lj) find. 

An interesting feature of Fig. El is the correlation be- 
tween the CS and SO abundances. Here the normalization 
to the CO abundance also serves as a valuable test: if for 
some reason the absolute density scale had been system- 
atically overestimated for the most massive envelopes and 
underestimated for the least massive envelopes, a false 
trend of abundances with mass could result and trends 
between abundances such as those seen in Fig. El should 
arise. In this case, however, normalization by a "standard" 
abundance should take out such an effect, but as illus- 
trated in Fig.^Jthis is not the case. The relation between 
CS and SO therefore seems to be real. 

Interestingly, the CS/SO abundance ratio has previ- 
ously also been suggested to trace evolutionary effects re- 
lated to cloud conditions and evolution, e.g., variation of 
the initial C/O ratio, density effects, the te mporal evolu- 
tion of a given core or imp ortance of X-rays (|Bergin et alJ 
Il997t iNilsson et all 120001) . It is found through time de- 
pendent modeling of the chemistry that the CS/SO ratio 
increases throughout the evolution of a molecular cloud 
starting from an atomic carbon-rich phase, but stabilizes 
at late times at a level dependent on the initial C/O ratio. 
As illustrated in Fig. El the relationship between the CS 
and SO abundances is clearly non-linear, implying that 
one or more of these effects may play a role in determin- 
ing the relative abundances of these two molecules. The 
CS/SO ratio varie s from sa 0.2 to 4, i n good agreement 
with the results of INilsson et al.l l|200f)() who analyzed CS 
and SO abundances from a sample of 19 molecular clouds. 

SO2 is detected toward only a few sources in the sam- 
ple. Typically, the upper limit to the SO2 abundance is 
found t o be a fewxlO" 10 in t his study. The same was 
seen by IBuckle fc Fuller! l)2003|) who only detected S0 2 
emission toward 30% of their sources, i.p ., sources in the 
Serpens region. In fact, IB uckle fc Fuller! did not detect 
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Fig. 13. CS vs. SO abundance. The dashed line indicates 
a linear relation between the CS and SO abundances, the 
solid line is the best-fit correlation. In the lower panel 
the abundances have been normalized to a CO abundance 
of 10~ 4 , mimicking the assumption in IBuckle fc Fuller! 
l|2003|) . Symbols are defined in Fig. ITU 



SO2 for any of the four sources also in our sample. For 
these sources up per limits based on the observations of 
IBuckle fc Fuller! are a fewxlO -11 . 

The non-dete c tions can also be compared to the results 
of ISchoier et"afl l|2002|) for IRAS 16293-2422, for which 
abundance jumps, either due to th ermal evapora tion or 
outflow-induced shocks, were found. ISchoier et al.l argued 
for an SO2 abundance jump from 2xl0 -10 in the outer 
envelope to lxlO -7 in the inner envelope. The SO2 lines 
in this study are in fact expected to probe the outer region 
of the envelope and the derived upper limit to the abun- 
dances do seem to indicate that the abundances found 
for IRAS 16293-2422 are higher than those found here. 
It is interesting to note that SO2 is only detected to- 
ward regions with high outflow activity (i.p., the objects 
in NGC 1333 and VLA1623) and that the SO 8 7 -7 6 was 
found to be very broad (and only detected) toward these 
objects with widths of ~ 5 — 10 km s _1 (FWHM) contrast- 
ing the other observed lines. These objects also show the 
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highest SO abundances. Together with the strong SO and 
SO2 emission toward the Serpens sources which are also 
related to strong outflows, this suggests an enhancement 
of sulfur-bearing species in the inner envelopes due to out- 
flows. Large enhancements of the sulfur-species (together 
with CH3OH and SiO) are observed in outflows where 
these c an be studied well separated from their driving pro- 
tostar l|Bachiller fc Perez GutierreJl997trjOrgensen et al.1 
20033)- A deep systematic study of the line emission from 
these and other sulfur species (e.g., H2S, HCS + , H2CS) to- 
ward a large sample of objects will shed more light on this 
question and thus provide better insight into the sulfur- 
chemistry in low-mass protostars. 

4.3. HCO+ and N 2 H + 

HCO + is of great importance in chemical models of pro- 
tostellar environments as it is the primary molecular ion 
and thus r egulator of the elect ron density/ionization struc- 
ture fe.g.. lCaselli et alJl2002() an d the most important de - 
stroyer of other molecules (e.g., iBergin fc Langerlll997l) . 

As shown in the upper panel of Fig. El the derived 
HCO + abundances sh ow an e volution with mass similar 
to that found for CO (|Paper J) . This is even more clearly 
illustrated in the lower panel of Fig.EJ where a tight cor- 
relation between CO and HCO + abundances is seen. In 
fact, the CO and HCO + abundances are linearly depen- 
dent with 

[HCO+] = 7.4 x 10~ 5 x [CO] 

or put differently: a "standard" undepleted CO abundance 
of 10 -4 corresponds to an HCO + abundance of 7.4 x 10~ 9 . 

It is found that N2H 4 " marks a clear contrast to HCO + : 
as shown in Fig. El the N2H" 1 " abundance decreases with 
increasing CO abundance. High angular res olution inter- 
ferometer maps of protostell ar regions (e.g., IBergin et al] 
l2OOlllj0rgensen et alJl2003bh find that cores with low CO 
abundances show up stronger when mapped in N2H + . 

Both trends can be understood when considering the 
chemical network in more detail taking the depletion of 
CO into account. For both HCO + and N2H + the primary 
formation routes are through reactions with H^~, i.e.: 

CO^HCO++H 2 (6) 

N 2 ->N 2 H++H 2 (7) 

For standard CO abundances ([CO/H 2 ] - 10" 4 ) eq. © 
is the dominant removal mechanism for Hjj~, but as CO 
freezes out this reaction drops in importance and eq. Q 
becomes more important for the removal of Hjj~ . The main 
destruction mechanism for HCO + is dissociative recom- 
bination, which is also the case for N2H" 1 " when CO is 
depleted. However, as CO returns to the gas-phase, de- 
struction of N2H" 1 " through reactions with CO: 



H 



H 



o 



10 











10~ 8 




:0 


□ 




♦ 

* ♦ 


_ 


_Q 

10 


=• 




♦ 

■ 

* 




10- 10 


< 


> 






0.1 


Mass / 


1.0 

' M 




10- 8 


r 


♦ 

O 

♦ 

♦ y 




O 


- 


10 J 


- 

♦ 


♦ ♦ 

O 






10- 10 










10~ 


5 


10 1 





N 2 H 4 



CO -> HCO^ 



N, 



(8) 



CO 



Fig. 14. HCO + abundance vs. mass (upper panel) and vs. 
CO abundance (lower panel). In the lower panel has the 
linear correlation between the HCO + and CO abundances 
been overplotted. Symbols defined as in Fig. IT21 

becomes the dominant removal mechanism for N 2 H + . In 
Appendix ^ we consider the chemical network for Hjj" , 
HCO + , and N 2 H + in detail. The main conclusions are that 
a linear increase of the HCO + abundance with CO abun- 
dance is expected when CO is depleted. For higher CO 
abundances, however, the HCO + abundance does not de- 
pend on [CO] since a balance between formation through 
eq. (JBJ and destruction through dissociative recombina- 
tion exists. In contrast the N 2 H + abundance is high when 
CO is depleted but declines rapidly as ([CO]) -2 with 
the increasing CO abundance as H^" is removed (forming 
HCO+) and N 2 H+ is destroyed through eq. I©. 

To further illustrate these points the upper panel of 
Fig. El shows the chemical network for low (depleted) 
and standard CO abundances. The lower panel shows 
the N 2 H + , HCO + and H^ abundances as functions of 
CO abundance calculated in a cell with density n(H2) = 
1 x 10 6 cm -3 and temperature T = 20 K at 10 4 years using 
the chemical code of S.D. Doty and adopting the chemistry 
used in the detailed chem ical modeling of the envelope 
around IRAS 16293-2422 ijDotv et al.ll2004h . The figure 
clearly shows the linear relationship between the CO and 
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HCO + abundances for CO values lower than w 2xl0 -5 
and likewise the rapid decline of N 2 H + for higher CO 
abundances. The absolute values of the abundances and 
the exact CO abundance dividing between the "low" and 
"standard" [CO] regions is regulated by the exact details 
of the chemistry (e.g., the initial N2 abundance) and the 
cosmic ray ionization rate, but the overall trends remain 
the same. Thus trends of a linear increase of HCO + abun- 
dance with increasing CO abundance can be understood 
in a limit where CO is depleted and Eq. © is no longer 
the dominant removal mechanism for H3 . 

4.4. HCN, HNC and CN 

HCN is the molecule with the most striking lack of cor- 
relation with mass or CO abundances, as can be seen 
in Fig. El Chemically HCN and its geometrical iso- 
mer, HNC, are naturally thought to be closely related 
and [HNC] /[HCN] ratios of unity or slight ly higher are 
typically observed toward molecular clouds (jHirota et alJ 
ll998HDickens et a.lJl20ort . 

For both HCN and HNC it is found that the 1-0 
lines trace material with higher abundances - or addi- 



Fig. 16. Upper panel: the chemical networks for low CO 
abundances (i.e., depletion) and standard CO abundance 
([CO]~ 1CP 4 ). The dominant reactions are indicated by 
solid arrows, secondary reactions by dashed arrows. Where 
dissociative recombination is the main destruction for a 
molecule (i.e., N2H + or HCO + ) this has been indicated 
by a dotted arrow. Lower panel: the electron, N2H+, H^", 
and HCO + abundances as functions of CO abundance in 



a cell with density n(H.2) 
T = 20 K. 



1 x 10 6 cm 3 and temperature 



tional material outside what can be described by the 
single power-law density models. As seen in Fig. 1171 the 
[H CN/H 13 CN] ratio is sig nificantly lower than 70 quoted 
by IWilson fc E.oodl (^994) . One exception is the case of 
N1333-I4B which has a high estimated HCN abundance, 
possibly related to confusion with the outflow. This ratio 
is, however, not correlated with mass, as would be ex- 
pected in case of an error in the opacity treatment of the 
lines. The explanation is more likely that the H 13 CN abun- 
dances are heavily biased toward determinations based 
on the low J lines observed with the Onsala telescope 
since the higher J lines are only detected toward a small 
fraction of the sources. Since the abundances derived on 
the basis of the isotopic H 13 CN thereby probe the outer- 
most, less depleted regions this should lower the estimated 
[HCN/H 13 CN] ratios. 
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A higher degree of CO depletion could be expected 
to lead to a removal of gas-phase carbon and oxygen and 
thereby a decline of the [HNC]/[HCN] and [CN]/[HCN] 
ratios. On the other hand it is found that neither the 
[CN]/[HCN] nor the [HNC]/[HCN] ratio correlate with the 
degree of CO depletion. Another option is destruction of 
HNC at higher temperatures through neutral-neutral re- 
actions. This would be in agreement with the result that 
the Orion molecular clouds have significantly lo wer HNC 
abundances relative to H CN llSchilke et al . 1992) than the 
dark clouds surveyed bv lHirota et alJ l)l998|) . 

Fig. E| illustrates the close correlation between the 
HNC and CN abundances also indicated by the correla- 
tion coefficients (Tabled and Fig. El- HNC and CN are 
expected to be related, with HCNH+ as an intermediate 
product, through the reactions: 



HNC 



H 3 + 



HCNH+ + H 2 



HCNH+ + e~ -> CN + H 2 



(9) 
(10) 



These reactions are according to the UMIST database 
(iLe Teuff et alJ l200(t the dominant formation and re- 
moval mechanisms for the three species at 20 K and 
lxlO 6 cm -3 . The main formation mehcanism for HCN at 
this temperature and density is also through dissociative 
recombination for HCNH+ but this is secondary compared 
to the formation of CN, which could explain the weaker 
correlation between HCN and the other nitrogen-bearing 
species. 

4.5. HQN 

The HC3N abundance has been suggested to be an indi- 
cator ofJ^MHgmEoraleyolution^ 

fg-g.jHira.ha.ra et alJIl 9Q2tlRuffle et alll QQTtlCaselli et alJ 
199q) in dark clouds and pre-stellar cores. The HC3N 
abundance peaks early in the evolution of dark clouds 
when a substantial amount of carbon is in atomic form 
in the gas-phase, but also increases with increasing deple- 
tion (i.e., potentially at "later" stages). Depletion tends 
to remove atomic oxygen from the gas-phase, which oth- 
erwise has a tendency to destroy ions necessary for the 
formation of species such as HC3N. Fig. 1201 compares the 
HC3N abundance with the CO abundance and the CS/SO 
abundance ratio. As can be seen, [HC3N] is not particu- 
larly higher in objects with a larger degree of CO depletion 
- except for the pre-stellar cores when these are considered 
separately (see Sect. 14.7(1 . 

For the protostars in our sample, however, Fig. 1201 
shows that the HC3N abundance is related to the 
[CS]/[SO] ratio - with lower ratios of the two sulfur- 
bearing molecules corresponding to lower HC3N abun- 
dances. This can be understood in a scenario where the 
HC3N abundance is indeed a tracer of atomic carbon, since 
the CS/SO ratio would likewise be increased by higher 
a mounts of atomic ca rbon, as suggested by the models 
of iBergin et al The question is then whether this 

should be taken as an indicator of chemical "youth". As 
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Fig. 17. The ratio of the CS and C 34 S abundances plot- 
ted vs. HCN and H 13 CN ratio. The big cross mark the 
predictions from the standard isotopic ratio of 12 C: 13 C of 
70 and 32 S: 34 S of 22. Symbols as in Fig.EJ 

can be seen in Fig. [20| the dynamically "older" class I ob- 
jects have higher HC3N abundances and [CS]/[SO] ratios, 
which apparently would contradict this suggestion. 

An alternative explanation could be that the amount 
of atomic carbon is enhanced by the impact of UV ra- 
diation from the outside due to the interstellar radiation 
field. We can, however, argue that this is not the case 
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from the CN line observations. As noted above the HC3N 
and CN abundances are found to be interlinked, which is 
not difficult to understand if one considers the dominant 
formation and destruction mechanisms for HC3N in gas 
where the degree of CO depletion is low: 

CN + C 2 H 2 -> HC 3 N + H 
C+ + HC 3 N -► C 3 H+ + CN 

If the trend seen between the HC3N and the CS/SO ra- 
tio is related to the radiation field, the correlation should 
be more clear for abundances constrained by the CN 1-0 
lines which probe the outermost region of the envelope. As 
illustrated in Fig. 1211 however, the trend is much stronger 
for the CN abundances constrained by the CN 3-2 lines 
that have a significantly higher critical density and thus 
probe the inner, more shielded region of the envelope. UV 
radiation from the central star could increase the amount 
of atomic carbon in the inner envelope and would lead to 
trends similar to that observed, but whether UV radia- 
tion can penetrate far enough out into the envelope to be 



important is still not clear. A more detailed study of the 
molecular species in the chemical network for the nitrogen- 
bearing species together with more detailed modeling in- 
cluding the radial variation of the molecular abundances 
is needed to fully address these questions. 



4.6. Deuterium fractionation 

The deuterium fractionation of HCO + is seen from plots 
of the [DCO+]/[HCO+] ratio in Fig.E2 The fact that the 
DCO + emission predominantly originates in the cold outer 
gas is evidenced by the narrow line widths for all sources: 
the turbulent broadening required to model the DCO + 
line widths is only 0.3-0.5 km s" 1 . The prestellar cores 
clearly show the highest [DCO+l /PHCO+1 ratio of ~ 5% in 
agreement with findings by, e.g.. lCaselli et alJ l)2002|) . The 
class sources show [DCO + ]/[HCO + ] ratios ranging from 
0.004 to 0.05. DCO + is not detected for the class I sources 
corresponding to the upper limits on the [DCO + ]/[HCO + ] 
ratio of ~ 0.001. As can be seen from the lower panel 
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Fig. 19. [CN] vs. [HCN] (upper panel) and vs. [HNC] 
(lower panel). Symbols as in Fig. 1121 



of Fig. the [DCO+]/[HCO+] ratio does seem to be 
correlated with the degree of CO depletion. 

The typical deuterium abundance ratios are in general 
significantly higher than the "cosmic" D/H ratio of 10 -5 . 
Both gas-phase reactions and grain-surface reactions have 
been invoked to describe the deuterium fractionation at 
low temperatures in pr e- and protostellar enviro nments. 
In gas-phase models bv lRoberts fc Millar! l|2000bh such a 
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Fig. 20. [HC3N] vs. mass and [CO] (upper panel) and vs. 
[CS]/[SO] ratio (lower panel). Symbols as in Fig. 1721 

trend is indeed expected. In pure gas-phase models the 
primary mechanism for driving the fractionation of HCO + 
is the small zero-point energy in the reaction: 



HT + HD *=> H?D 



H 2 



(11) 



that predominantly drives the D into H2D + relative to H3 , 
and which subsequently reacts with CO to form DCO + . 
Depletion of CO causes (jlljl to be the dominant mecha- 
nism for removal of and since DCO + is produced sub- 
sequently, the formation of HCO + (eq. © in Sect. l4.3|l will 
be less pro ductive compared to the deute rated versions. 
Mode ls by iRoberts fc Millarl l|2000a|) and iRoberts et all 
(200!§) show that a high degree of depletion may be re- 
quired in order to produce formation of the doubly and 
triply deuterated species observed in protostellar environ- 
ments llCeccarelli et alJll998t |Lis et alJl2002t IFarise et alJ 



120021: Ivan der Tak et al.ll2002h . Other studies indicate an 
increase of the deuteration with CO depletion, e.g., the 
survey of the deuteration of Do CO f or a sample of pre- 
stellar cores bv lBacmann et a~ 

Figure E3| compares the [DCN]/[HCN] ratio with en- 
velope mass and [DCO + ]/[HCO + ] ratio. It appears that 
the [DCN]/[HCN] and [DCO+]/[HCO+] ratios are not cor- 
related. In fact the absence of DCN in the pre-stellar 
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Fig. 21. CS/SO abundance ratio vs. abundance of CN 
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probing the outer and inner regions of the envelope, re- 
spectively. Symbols as in Fig. 



cores is striking considering their strong DCO + emission. 
This could imply that the [DCN]/[HCN] ratio is higher 
for the warmer envelopes. This is in agreement with gas- 
phase deuteration of HCN, which may occur at slightly 
higher temperatures than that of HCO + : in particular 
deuteration through CH^" + HD — * CH2D+ may be more 
important fo r temperatures higher than s» 30 K (e.g., 
iTurnerlEooH) . Alternatively, the [DCN]/[HCN] ratio may 
have been established earlier in the protostellar evolution, 
frozen out onto the dust grains and released back at higher 
temperatures than is the case for the [DCO + ]/[HCO + ] ra- 
tio. 

4. 7. The pre-stellar cores 

The two pre-stellar cores in our sample , L1689B and 
LI 544, were also studied b vlLee et al. I ll2003h together with 
an additional core, L1512. lLee et alJ found a high degree of 
depletion of CO and HCO + for L1544 but close to "stan- 
dard" abundances for these molecules in L1689B. This is 
consistent with the results in this paper. iLee et alJ also 



Fig. 22. [DCO + ]/[HCO + ] ratio vs. mass (upper panel) vs. 
and [CO] (lower panel). Symbols as in Fig. 1121 



observed N2H+ 1-0 and found it to be weak in L1689B 
compared to L1544. This is in good agreement with our 
results which show N2H + to be an order of magnitude 
more abundant in L1544 than in L1689B, strengthening 
the N2H + and CO anti-correlation discussed in Sect. 14.31 
Of the other molecules in this paper, the sulfur-bearing 
species (CS and SO) are also found to have lower abun- 
dances in L1544 indicating a higher degree of overall de- 
pletion. The nitrogen-bearing species in contrast show an 
opposite trend with high abundances in L1544. In par- 
ticular, HC3N is close to a factor of 50 higher in L1544 
than in L1689B, supporting the suggestion that HC3N 
traces the degree of depletion in the p re-stellar stages (e.g. 
iRuffle et al.lll997tlCaselli et al.lll998l) . Also the higher de- 
gree of HCO + deuteration in L1544 than in L1689B is 
consistent with the higher degree of depletion L1544 than 
in L1689B (ref. the discussion in Sect. 14. 6[) . 

4.8. Comparison to other star-forming regions 

In Fig. ^1] the average abundances for the class 
and I objects and pre-stellar cores are compared to 
the abundances of other star-forming regions, the outer 
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upper limit found for the objects in this study. It does 
not, however, stand markedly out considering the scat- 
ter in abundances within the larger group of class ob- 
jects. It therefore seems that IRAS 16293-2422, despite 
possibly being affec ted by outflows on smaller scales (e.g., 



Schoier^t^alji 2004) a nd having a "hot in ner region" (e.g., 



Ceccarelli et al . 1998; Schoi er et al . 2002) has a cold outer 



envelope that is similar to that of the other class objects 
in terms of the overall abundances. Otherwise the most 
striking feature of Fig. Ellis the significantly higher HNC 
abundance (two orders of magnitude) in L134N compared 
to the other sources and molecules. The high- and the low- 
mass YSOs differ slightly with SO and HCN abundances 
higher by up to a factor 5 and HC3N lower by a factor 
5-10 in the high-mass YSOs. 
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envelope around th e class object IRAS 16293-2422 
dSchoier et, al.ll2002l) . the average abundance for the 3 



(Helmich & van Dishoeckl 


Il997l) and the "( 


the dark cloud L134N 


Dickens et alJ 


2000) 



As mentioned in the introduction the class object 
IRAS 16293-2422 is the most studied low-mass proto- 
star in terms of the chemistry of its protostellar en- 



velope (see, e.g. . IB lake etal.1 Il994l Ivan Dishoeck et all 
19951 ICeccarelli et al.lll998l l2000albl: ISchoier et alJl2002t 
P arise et alJl2002t ICazaux et alJl2003h . because of its rich 
spectrum and its warm inner region where ices have 
evaporated. This naturally raises the question whether 
IRAS 16293-2422 is indeed a typical class object: is the 
richness of its spectrum simply caused by it being the 
closest object with the most massive envelope, or is it 
caused by other effects, such as the interaction of its out- 
flow with the nearby envelope? As it can be seen from 
Table [T^l and Fig. 1121231 IRAS 16293-2422 has a fairly 
standard set of outer envelope abundances for CO, CS 
and HCN. On the other hand it shows lower abundances 
(factors 4-20) of especially HNC, CN and N 2 H+ and high 
abundances of SO and SO2 compared to the typically 



5. Conclusion 

The molecular inventories for the envelopes around a 
sample of low-mass protostars have been established. 
Using models for the one dimensional physic al structure 
of the envelopes from Ijprgensen et alJ l|2002l) (Paper I), 
the abundances of a range of molecular species are con- 
strained through Monte Carlo line radiative transfer mod- 
eling of single-dish submillimcter and millimeter observa- 
tions. The main conclusions are: 

1. For most sources and molecules, the high excitation 
0.6-1.4 mm lines are well represented by a constant 
fractional abundance. The low excitation 3 mm lines of 
HCO + and the nitrogen-bearing species, however, are 
signifcantly underestimated by such models, similar to 
the trend seen for CO in Paper I. Varying freeze-out 
timescales in the regions of the envelopes correspond- 
ing to different densities can explain this. All lines for 
these species can be accounted for with a "drop" profile 
where the envelopes have "standard" molecular abun- 
dances in the outermost low density region and in the 
innermost high temperature region, but drops in an 
intermediate zone where the molecule can freeze-out. 

2. Effects of the envelope velocity field are tested and con- 
firmed to be insignificant when determining the abun- 
dances from optically thin species. An upper limit to 
the magnitude of the velocity field and thereby the 
mass accretion rate can be inferred from comparison 
to the widths of those lines. 

3. An empirical chemical network is constructed through 
calculations of correlation coefficients for pairs of abun- 
dances. Strong correlations are found between, e.g., 
HCO + and CO, CS and SO and between the nitrogen- 
bearing species, HNC, CN and HC 3 N. 

4. A linear relationship between the abundances of CO 
and HCO + is seen, whereas the N2H" 1 " abundances ap- 
pear lower in the objects with high CO abundances. 
This can be understood from analytic considerations 
of the chemical network of HCO + and N2H + , where 
both are formed through reactions involving H^. A 
linear relationship between [CO] and [HCO + ] exist in 
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the limit where N2 rather than CO removes H3 due to 
CO freeze-out. Likewise the N 2 H+ abundances rapidly 
decline as soon as the CO abundances increase, be- 
cause and N 2 H + are destroyed through reactions 
with CO. 

5. The CS/SO abundance ratio is found to be correlated 
with the abundances of CN and HC3N, which may 
reflect the amount of atomic carbon in the gas-phase. 

6. High levels of deuterium fractionation of HCO + 
([DCO+]/[HCO+] = 1-5 %) are found in the pre- 
stellar cores and in the most massive/coldest envelopes 
around the class objects. The [DCN]/[HCN] ratio is 
not correlated with the [DCO + ]/[HCO + ] ratio, possi- 
bly indicating differences in the gas-phase deuteration 
of HCO+ and HCN. 

7. Among the two pre-stellar cores in our sample, L1544 
shows a higher degree of depletion of CO and the 
sulfur-bearing species, resulting in a chemistry with 
higher abundances of N 2 H + , HC3N and a higher 
degree of HCO + deuteration than what is seen in 
L1689B. 

8. The outer abundances of the p reviously studied cla ss 
object IRAS 16293-2422 fe.g.. lSchoier et al.ll2002h are 
found to be in agreement with the average abundances 
in the studied sample of class objects, indicating that 
the chemistry of the outer cold envelope of this partic- 
ular object is similar to that of other class objects. 

The results of this paper show the potential of multi- 
transition line observations to understand radial vari- 
ations and find empirical correlations between various 
species elucidating the chemistry in star-forming envi- 
ronments. Future observational studies should focus on 
deeper systematic searches, e.g., in order to fully under- 
stand the sulfur chemistry and to constrain abundances 
of species connecting the nitrogen-bearing species to the 
overall chemical network. Also high-resolution observa- 
tions with, e.g., the SMA, CARMA and ALMA can be 
used to test the p redictions of radial variations; as dis- 
cusse d elsewhere l|j0rgensen et all l2003bl ISchoier et all 
12004^1 the models presented here serve as a valuable start- 
ing point for interpreting such data. This study has not 
addressed the chemistry of the innermost hot core region 
where complex organic molecules may be formed through 
evaporation of ices and high-temperature gas-phase chem- 
istry. Since these regions are small (typically < 1" for the 
objects considered here), high angular resolution observa- 
tions of high-excitation lines are needed to probe them. 
Finally, the results presented here will form the basis for 
more detailed modeling of the envelope chemistry. 
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whereas for "standard" CO abundances of 1x10 4 the 
main removal mechanism becomes: 



N 2 H+ + CO -> N 2 + HCO H 



(A.5) 



In the low [CO] limit the main removal mechanism for 
is by eq. (|A. 31) rather than by eq. (|A.lfl . so in this limit 
we write for the formation/destruction balance for 



For H+ the main formation route is through cosmic ray 
ionization of H2 : 

H 2 + CR -> H+ + e" 



<>h 2 = fc N2H+ n H+ n N2 



(A.14) 



replacing eq. (|A.8|1 . From eq. (|A.7|I and (IA. 14|) we can then 
eliminate nnco+ : 



«HCO + 



Cm 2 



H 2 



H 



H 



(A.6) 



3 is through eq. (|A.3|I 
when the CO abundance is low and through eq. (|A.1(I 



The main removal mechanism for H3 is through eq 



"-HCO+ = n CO- , 

kco n e % 2 h+ nn 2 

or in terms of the abundances: 

&HCO+ C/™H 2 



[HCO^ 



«CO «N 2 H+ 



[e] [N 2 



[CO] 



(A.15) 



(A.16) 



when the CO abundance is close to the standard value 
~ 10 4 . 

We introduce the rate coefficients for each reaction 
with fcjjco+ j ^COi an d fci\i 2 H+ being the rate coefficients 
for eq. (|A.1(1 . (|A.2|I . and (|A.3|I . respectively, and /cn 2 and 
fc^ 2 being the rate coefficients for eq. (|A.4(1 and l|A.5|> . 
respectively. 

We can equate the formation and destruction rates of and again eliminating n H + , rewriting the expression in 



Furthermore, in this limit the destruction of N 2 H + is by 
dissociative recombination (eq. (|A.4|I . so for N 2 H + we can 
write: 



&n 2 h+ »%+ n N 2 = fe N 2 1c% 2 H+ 



(A.17) 



HCO+ as: 

fcHco+ n u+ n co = k co «hco+ n e 



terms of abundances, we find: 
C/™H 2 



(A-7) [N2H+] = 



where n HC o+ , n H + , nco and n c are the densities of HCO + , 
, CO and electrons, respectively. 



/cn 2 [e] 



(A.18) 



The main points to be drawn from these equations 
(i.p., eq. l|A~TUl) . l|A~llij) . (|A~TB|l and ||A~T%1) ) are that the 



When the CO abundance is - 10 4 the main removal HC0 + abundance does not depend on the CO abundance 
mechanism for H+ is eq. (53) so equating the formation in the « high » ^ whereas the N2 h+ abundance de- 



and destruction rate for H3 gives 

Oh 2 = fcHCO+ n H+ n CO 



(A. 



clines quickly as ([CO]) -2 . In the low [CO] limit, however, 
the HCO + abundance increases linearly with increasing 
CO abundance whereas the N 2 H + abundance is constant. 



where £ is the cosmic ray ionization rate, ?ih 2 the H 2 den- 
sity and it is assumed that each molecular ion pro- 
duced in eq. (|A.6|I immediately reacts with H 2 to form 



(A.9) 



. From eq. (S3 and JSjj) we can then write: 
Cn H2 = fcco n H co+ n e 



which dividing by n\ , introducing the abundances [X] = 
nx/nn 2 and isolating [HCO + ] gives: 



[HCO^ 



C/"H 2 

kco N 



(A.10) 



We can similarly equate the formation and destruction 
rates for N 2 H+ (eq. (|A.3j) and HA.5|l ): 

fc N 2 H+" H + n N 2 = k' N2 n N2li +n C o ( A -H) 
which together with eq. (|A.8|) gives (eliminating n H +): 

Oh 2 



^N 2 H+ n N 2 



k' 



N 2 n N 2 H+ "CO 



(A.12) 



&HCO + "CO 

which can be rewritten in terms of the N 2 H + abundance: 
[N 2 ] fc N2 H+ 



[N 2 ET 



[CO] 2 k' fc H co- 



■ C/«H 2 



(A.13) 



